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Abstract 

A Delaunay cell decomposition of a surface with constant curvature gives rise to a 
circle pattern, consisting of the circles which are circumscribed to the facets. We treat 
the problem whether there exists a Delaunay cell decomposition for a given (topological) 
cell decomposition and given intersection angles of the circles, whether it is unique and 
how it may be constructed. Somewhat more generally, we allow cone-like singularities 
in the centers and intersection points of the circles. We prove existence and uniqueness 
theorems for the solution of the circle pattern problem using a variational principle. The 
functionals (one for the euclidean, one for the hyperbolic case) are convex functions of 
the radii of the circles. The critical points correspond to solutions of the circle pattern 
problem. The analogous functional for the spherical case is not convex, hence this case is 
treated by stereographic projection to the plane. From the existence and uniqueness of 
circle patterns in the sphere, we derive a strengthened version of Steinitz' theorem on the 
geometric realizability of abstract polyhedra. 

We derive the variational principles of Colin de Verdiere, Bragger, and Rivin for circle 
packings and circle patterns from our variational principles. In the case of Bragger's and 
Rivin's functionals, this requires a Legendre transformation of our euclidean functional. 
The respective Legendre transformations of the hyperbolic and spherical functionals lead 
to new variational principles. The variables of the transformed functionals are certain 
angles instead of radii. The transformed functionals may be interpreted geometrically as 
volumes of certain three-dimensional polyhedra in hyperbolic space. Leibon's functional 
for hyperbolic circle patterns cannot be derived from our functionals. But we construct 
yet another functional from which both Leibon's and our functionals can be derived. By 
applying the inverse Legendre transformation to Leibon's functional, we obtain a new 
variational principle for hyperbolic circle patterns. 

We present Java software to compute and visualize circle patterns. 

Zusammenfassung 

Eine Delaunay-Zellzerlegung einer Flache konstanter Kriimmung liefert ein Kreismus- 
ter, welches aus den Kreisen besteht, die den Facetten umschrieben sind. Wir betrachten 
das Problem, ob es fiir eine vorgegebene (topologische) Zellzerlegung und vorgegebene 
Schnittwinkel zwischen den Kreisen eine entsprechende Delaunay-Zellzerlegung gibt, ob 
sie eindeutig ist, und wie sie zu konstruieren ist. Etwas allgemeiner lassen wir auch ke- 
gelartige Singularitaten in den Mittel- und Schnittpunkten der Kreise zu. Wir beweisen 
Existenz- und Eindeutigkeitssatze fiir die Losung des Kreismusterproblems mit Hilfe von 
Variationsprinzipien. Die Funktionale (eins fiir den euklidischen, eins fiir den hyperboli- 
schen Fall) sind konvexe Funktionen der Radien der Kreise. Kritische Punkte entsprechen 
Losungen des Kreismusterproblems. Das analoge Funktional fiir den spharischen Fall ist 
nicht konvex, deshalb wird dieser Fall durch stereographische Projektion in die Ebene 
erledigt. Aus der Existenz und Eindeutigkeit von Kreismustern in der Sphare folgern wir 
eine verscharfte Version des Satzes von Steinitz iiber die geometrische Realisierbarkeit von 
abstrakten Polyedern. 

Wir leiten die Variationsprinzipien von Colin de Verdiere, Bragger und Rivin fiir 
Kreispackungen bzw. Kreismuster aus unseren Variationsprinzipien ab. Im Fall der Funk- 
tionale von Bragger und Rivin erfordert dies eine Legendretransformation unseres eukli- 
dischen Funktionals. Entsprechende Legendretransformationen des hyperbolischen und 
des spharischen Funktionals liefern neue Variationsprinzipien. Die Variablen der transfor- 
mierten Funktionale sind nicht Radien, sondern bestimmte Winkel. Die transformierten 
Funktionale besitzen eine geometrische Interpretation als Volumen von bestimmten drei- 
dimensionalen Polyedern im hyperbolischen Raum. Leibons Funktional fiir hyperbolische 
Kreismuster lasst sich nicht aus unseren Funktionalen herleiten. Wir konstruieren jedoch 
ein weiteres Funktional, aus dem sowohl Leibon's als auch unser Funktional hergeleitet 
werden kann. Durch die umgekehrte Legendretransformation von Leibons Funktional er- 
halten wir ein neues Variationsprinzip fiir hyperbolische Kreismuster. 

Wir prasentieren Java Software zur Berechnung und Visualisierung von Kreismustern. 
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CHAPTER 1 



Introduction 

1.1. Existence and uniqueness theorems 

A circle packing is a configuration of circular discs in a surface sucli that tlie 
discs may touch but not overlap. We consider only circle packings in surfaces 
of constant curvature. Connecting the centers of touching discs by geodesies as 
in figure 11.11 one obtains an embedded graph, the adjacency graph of the packing. 
Consider the case when the adjacency graph triangulates the surface as in figure [TTTl 
(left). The following theorem of Koebe j26j answers the question: Given an abstract 
triangulation of the sphere, does there exist a circle packing whose adjacency graph 
is a geometric realization of the abstract triangulation? 

Theorem 1.1 (Koebe). For every abstract triangulation of the sphere there is 
a circle packing whose adjacency graph is a geometric realization of the triangula- 
tion. The circle packing corresponding to a triangulation is unique up to Mobius 
transformations of the sphere. 

Now consider the case when the adjacency graph of a circle packing gives rise 
to a cell decomposition whose faces are not necessarily triangles, as in figure 11.11 
(right). To characterize the cell decompositions of the sphere which correspond to 
circle packings, we need the following definition. 

Definition. A cell complex is called regular if the characteristic maps, that 
map closed discs onto the closed cells, are homeomorphisms. A cell complex is 
called strongly regular if it is regular and the intersection of two closed cells is 
empty or equal to a single closed cell. 

Remark. Suppose a cell complex S is in fact a cell decomposition of a compact 
surface without boundary. One obtains the following conditions for E being regular 
and strongly regular. 

The cell decomposition E is regular if and only if the following conditions hold. 

(i) Each edge is incident with two vertices. (There are no loops.) 




Figure 1.1. Left: A circle packing corresponding to a triangulation. Right: A pair 
or orthogonally intersecting circle packings corresponding to a cell decomposition. 
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(ii) Each edge is incident with two faces. (There are no stalks.) 
(in) If a vertex v and a face / are incident, there are exactly two edges incident 
with both V and /. 

The cell decomposition S is strongly regular if and only if, in addition, the 
following conditions hold. 

(iv) No two edges are incident with the same two vertices. 

(v) No two edges are incident with the same two faces. 

(vi) If each of two faces is incident with each of two vertices, then there is an 
edge which is incident with both faces and both vertices. 

The above characterization implies: The cell decomposition E is (strongly) 
regular if and only if its Poincare-dual decomposition is (strongly) regular. 

A cell decomposition of the sphere which arises from a circle packing is strongly 
regular. Conversely, Koebc's theorem implies that every strongly regular cell de- 
composition of the sphere comes from a circle packing. (Simply triangulate each 
face by adding an extra vertex inside and connecting it to the original vertices.) 
However, the corresponding packing is generally not unique up to Mobius trans- 
formations. The following theorem is a generalization of Koebe's theorem which 
retains uniqueness by requiring the existence of a second packing of orthogonally 
intersecting circles as in figure [lTI (right). 

Theorem 1.2. For every strongly regular cell decomposition of the sphere, there 
exists a pair of circle packings with the following properties: The adjacency graph 
of the first packing is a geometric realization of the given cell decomposition. The 
adjacency graph of the second packing is a geometric realization of the Poincare 
dual of the given cell decomposition. Therefore, to each edge there correspond four 
circles which touch in pairs. It is required that these pairs touch in the same point 
and intersect each other orthogonally. 

The pair of circle packings is unique up to Mobius transformations of the sphere. 

In the case of a circle packing corresponding to a triangulation as in Koebe's 
theorem, the second orthogonal packing always exists. Thus, theorem II. 21 is indeed 
a generalization of Koebe's theorem. The first published proof is probably due to 
Brightwell and Scheinerman |13| . They do not claim to give the first proof. The 
works of Thurston |45| and Schramm |40| (see theorem 11.41 below) indicate that 
the theorem was well established at the time. 

Associated with a polyhedron in E,"^ is a cell decomposition of the sphere rep- 
resenting its combinatorial type. We say that the polyhedron is a (geometric) 
realization of the cell decomposition. Steinitz' representation theorem for convex 
polyhcdra in R'^ says that a cell decomposition of the sphere represents the com- 
binatorial type of a convex polyhedron if and only if it is strongly regular |43| . 
|44j. Theorem 1 1.21 implies the following stronger representation theorem for convex 
polyhcdra in (see Ziegler |47j . theorem 4.13 on p. 118). 

Theorem 1.3. For every strongly regular cell decomposition of the sphere there 
is a convex polyhedron in M'^ which realizes it, such that the edges of the polyhe- 
dron are tangent to the unit sphere. Such a geometric realization is unique up to 
projective transformations which fix the sphere. 

Simultaneously, there is a polyhedron with edges tangent to the sphere which 
realizes the Poincare-dual cell decomposition such that corresponding edges of the 
two polyhcdra intersect each other orthogonally and touch the sphere in their point 
of intersection. 

Among the projectively equivalent polyhedral realizations, there is one and up 
to isometrics only one with the property that the barycenter of the points where the 
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Figure 1.2. A Dclaunay type circle pattern with Delaunay and Dirichlet cell de- 
compositions. 

edges touch the sphere is the center of the sphere. Every topological symmetry of 
the cell decomposition corresponds to an isometry of this polyhedral realization. 

The following much more general theorem is due to Schramm |40| . The proof 
is based on topological methods which are beyond the scope of this thesis. 

Theorem 1.4 (Schramm). Let P be a (3-dimensional) convex polyhedron, and 
let K C M'^ he a smooth strictly convex body. Then there exists a convex polyhedron 
Q C M"^ combinatorially equivalent to P whose edges are tangent to K . 

The two circle packings of theorem 11.21 form a pattern of orthogonally inter- 
secting circles. By way of a further generalization, one may consider circle patterns 
with circles intersecting at arbitrary angles. 

'Definition'. A circle pattern is a configuration of circles in a constant curva- 
ture surface which corresponds in some way to a cell decomposition of the surface. 
According as the constant curvature is positive, zero, or negative, we speak of 
spherical, euclidean, or hyperbolic circle patterns. 

To obtain a real definition, the correspondence between circle pattern and cell 
decomposition has to be specified. We will only be concerned with a special class 
of circle patterns which are connected to Delaunay cell decompositions. To be 
precise, we call them 'Delaunay type circle patterns', but we will usually refer to 
them simply as 'circle patterns'. Figure [L^ shows an example. 

Definition. A Delaunay decomposition of a constant curvature surface is a 
cellular decomposition such that the boundary of each face is a geodesic polygon 
which is inscribed in a circular disc, and these discs have no vertices in their interior. 
The Poincarc-dual decomposition of a Delaunay decomposition with the centers of 
the circles as vertices and geodesic edges is a Dirichlet decomposition (or Voronoi 
diagram) . A Delaunay type circle pattern is the circle pattern formed by the circles 
of a Delaunay decomposition. More generally, we allow the constant curvature 
surface to have cone-like singularities in the vertices and in the centers of the circles. 
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Figure 1.3. Interior intersection angle 9* and exterior intersection angle of two 
circles. 6 + 9* = tt. 

Figure [L^ shows a Delaunay decomposition (black vertices and solid lines), the 
dual Dirichlet decomposition (white vertices and dashed lines) and the correspond- 
ing circle pattern. The faces of the Delaunay decomposition correspond to circles. 
The vertices are intersection points of circles. 

Remark. In section ITTl we will give an alternative and slightly more general 
definition for Delaunay type circle patterns. 

A different class of circle patterns, Thurston type circle patterns, has been in- 
troduced by Thurston )45j . Here, the circles correspond to the faces of a cell 
decomposition, but the vertices do not correspond to intersection points. All ver- 
tices have degree 3. Two circles corresponding to faces which share a common edge 
intersect (or touch) with an interior intersection angle 9* (see figure n~5|) satisfying 
0<9* < 7r/2. 

Those Thurston type circle patterns in the sphere with the property that the 
sum of the three angles 9 around each vertex is less or equal to 27r correspond to 
hyperbolic polyhedra of finite volume with dihedral angles at most tt/2. Such poly- 
hedra have been studied by Andreev U, Pj. Delaunay type circle patterns in the 
sphere (without cone-like singularities) correspond to convex hyperbolic polyhedra 
with finite volume and all vertices on the infinite boundary of hyperbolic space. 

From a Delaunay type circle pattern, one may obtain the following data. (For 
the definition of interior and exterior intersection angle of two circles see figure [OH 
We will always denote the exterior intersection angle by 9 and the interior intersec- 
tion angle by 61*. Note that 9 + 9* ^n.) 

• A cell decomposition E of a 2-dimensional manifold. 

• For each edge e of S the exterior (or interior) intersection angle 9^ (or 91). 
It satisfies < 9^ < tt. 

• For each face / of E the cone angle <&/ > of the cone-like singularity 
at the center of the circle corresponding to /. (If there is no cone- like 
singularity at the center, then $/ = 27r.) 

Note that the cone angle at a vertex w of E is already determined by the 
intersection angles 9e. 

= E^e, (1.1) 
where the sum is taken over all edges around v. (See figure H"!^ ) If Qy = 2tt, there 
is no singularity at v. The curvature in a vertex v is 

= 2^-e„ (1.2) 

and the curvature in the center of a face / is 

Kf = 27r-^f. (1.3) 
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The following theorems assume that the surface is closed, and that there are no 
cone-like singularities. Only the last theorem II . 81 deals with the general case. 

Consider the following problem: Given such data as listed above, is there a 
corresponding circle pattern? If so, is it unique? The following theorem of Rivin |35| 
gives an answer for circle patterns in the sphere without cone-like singularities. 

Theorem 1.5 (Rivin, circle pattern version). Let T, he a strongly regular cell 
decomposition of the sphere and let an angle Of, with < 9^ < tt be given for every 
edge e ofE. Let E* be the Poincare dual decomposition of E, and for each edge e 
of Ti denote the dual edge of E* by e* . 

A Delaunay pattern corresponding to E with exterior intersection angles Og 
exists if and only if the conditions (i) and (ii) are satisfied. 

(i) If some edges e*, . . . , e* form the boundary of a face of E*, then 

(ii) If some edges e^^, . . . , e* form a closed path ofE* which is not the boundary 
of a face, then 

// it exists, the circle pattern is unique up to Mobius transformations of the 
sphere. 

The paper cited above deals with polyhedra in hyperbolic 3-space with vertices 
on the infinite boundary. The above theorem has the following, equivalent, form. 

Theorem 1.6 (Rivin, ideal hyperbolic polyhedra version). Let Y, be a strongly 
regular cell decomposition of the sphere and let an angle 9e with < 6^ < it be 
given for every edge e ofT,. Let E* be the Poincare dual decomposition ofY,, and 
for each edge eofYi denote the dual edge of E* by e* . 

A polyhedron in hyperbolic 3-space with vertices on the infinite boundary which 
realizes E and has exterior dihedral angles 6^ exists if and only if the conditions (i) 
and (ii) are satisfied. 

(i) If the edges , . . . , e* form the boundary of a face of T,* , then 

(ii) If the edges e|, . . . , e* form a closed path of E* which is not the boundary 
of a face, then 

^ > 2tt. 

If it exists, the polyhedron is unique up to an isometry of hyperbolic space. 

Theorem 11.71 below is a generalization of Rivin's theorem 11.51 to higher genus 
surfaces. Only the case of oriented surfaces has to be considered: Because of the 
uniqueness claim, non-orientable surfaces may be treated by applying the theorem 
to the orientable double cover. 

In the higher genus case, it is too restrictive to allow only strongly regular cell 
decompositions. For example, figure 11.41 shows a cell decomposition of the torus 
which is not even regular and the corresponding circle pattern with orthogonally 
intersecting circles. With appropriate cone-like singularities in the vertices and 
centers of circles, even faces whose boundary is a loop and vertices of degree one 
are possible. As a consequence of this tolerance regarding which cell decompositions 
are acceptable, theorems 1 1 . 71 and 11.81 are really only true for Delaunay type circle 
patterns in the sense of the slightly more general definition in section 12.11 

The condition of theorem 11.71 does not involve cycles in the Poincare dual 
decomposition like theorem ll.5l Instead, one has to consider cellular immersions of 
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Figure 1.4. Left: A cell decomposition of a torus (solid lines and black dots) and 
its Poincare dual (dashed lines and white dots). The top and bottom side, and the 
left and right side of the rectangle are identified. Right: The corresponding double 
periodic circle pattern with orthogonally intersecting circles. 

discs, whose interior is embedded, into the Poincare dual decomposition. The image 
of the boundary path of the disc is a closed path in which some edges might appear 
twice. For example, the shaded area in figure ^31 (l^ft) is the image of such a disc, 
where the heavy dashed lines are the image of its boundary. One edge appears 
twice in the image of the boundary path. The condition of theorem 11.71 is that 
the sum of over such a boundary, where the edges are counted with appropriate 
multiplicities, is at least 27r. 

Theorem 1.7. Let Yi he a cell decomposition of a closed compact oriented 
surface of genus g > 0. Suppose exterior intersection angles are prescribed by a 
function 9 S (0,7r)^ on the set E of edges. Then there exists a circle pattern 
corresponding to these data in a surface of constant curvature ( equal to if g = 1 
and equal to —1 if g > 1), if and only if the following condition is satisfied. 

Suppose A —f T,* is any cellular immersion of a cell decomposition of the closed 
disc A into the Poincare dual E* o/E which embeds the interior of A. Let ei, . . . , efe 
be the boundary edges of A, let be e^, . . . , e^. their images in S*, and let ei, . . . , Ck 
be their dual edges in S. (An edge of Y, may appear twice among the Cj.) Then 

k 

J20e,>27T, (1.4) 

where equality holds if and only if A has only one face. 

In the case they exists, the constant mean curvature surface and the circle 
pattern are unique up to similarity, if g = 1, and unique up to isometry, if g > 1. 

Schlenker j39j independently proves an existence and uniqueness result for hy- 
perbolic manifolds with polyhedral boundary, all vertices at infinity, and prescribed 
dihedral angles. Theorem 1 1 . 71 follows as a special case. Interestingly, to obtain the 
general result, Schlenker needs to first prove this special case separately. 

We deduce theorem fl . 7l from the following more technical, but also more general 
theorem 11.81 It is not assumed that 9 sums to 27r around each vertex. Hence 
there may be cone-like singularities in the vertices with cone angle 9 given by 
equation Hl.l|l . Also, cone-like singularities with prescribed cone angle are 
allowed in the centers of circles. Furthermore, it applies also to surfaces with 
boundary. For a boundary face /, the angle does not prescribe a cone angle, but 
the angle covered by the neighboring circles, as shown in figure [T31 (right). These 
angles on the boundary constitute Neumann boundary conditions. Alternatively, 
one might also prescribe the radii of the boundary circles. This would constitute 
Dirichlet boundary conditions. We consider only the Neumann problem. 
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Figure 1.5. Left: A cell decomposition of the disc. Right: A corresponding circle 
pattern with orthogonally intersecting circles. The shaded angles are the angles $ 
for boundary faces. Here, they are Stt/G for the three corner circles and 37r/2 for 
the other boundary circles. 



Theorem 1.8. Let be a cell decomposition of a compact oriented surface, 
with or without boundary. Suppose (interior) intersection angles are prescribed by 
a function 0* e (0,7r)^" on the set Eq of interior edges. Let $ £ (0,oo)^ be a 
function on the set F of faces, which prescribes, for interior faces, the cone angle 
and, for boundary faces, the Neumann boundary conditions. 

(i) A euclidean circle pattern corresponding to these data exists if and only if 
the following condition is satisfied: 

If F' Q F is any nonempty set of faces and E' C Eq is the set of all interior 
edges which are incident with any face f £ F' , then 



where equality holds if and only if F' ~ F. 

If it exists, the circle pattern is unique up to similarity. 

(ii) A hyperbolic circle pattern corresponding to these data exists, if and only if, 
in the above condition, strict inequality holds also in the case F' ^ F. If it exists, 
the circle pattern is unique up to hyperbolic isometry. 

Similar results were obtained by Bowditch 0, Garrett |19| . Rivin \S6\ . and 
Leibon |28| . Bowditch treats the euclidean case for closed surfaces without cone- 
like singularities in the centers of the circles. His proof is topological in nature: It 
hinges on the fact that a certain function (essentially the gradient of our functional 
S'cuc) is injective and proper. Rivin extends this result to surfaces with boundary. 
Also, he considers not only singular euclidean structures, but singular similarity 
structures. That is, he admits not only cone-like singularities (with rotational ho- 
lonomy) but singularities with dilatational and rotational holonomy. The proof uses 
his variational principle |34j . Leibon treats the hyperbolic case for closed surfaces 
and without cone-like singularities. The proof uses his variational principle (see 
section 1^21 • Garrett |19j obtains a similar theorem for euclidean and hyperbolic 
circle packings with prescribed cone-like singularities. He considers the Dirichlet 
boundary value problem (prescribed radii at the boundary). His proof uses the 
relaxation method developed by Thurston |45 |. 




(1.5) 



/GF' 
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1.2. The method of proof 

Chapter El contains our proofs of the theorems 11.21 11.31 11.51 11.61 11.71 and 11.81 

Here, we give a brief outhne of the proofs of theorems 11.81 11.71 and 11.51 which 
are the most involved. Most of the effort goes into proving the fundamental the- 
orem ^21 from which the others are deduced. We extend methods introduced by 
Colin de Verdiere for circle packings |16j . 

First, the geometric problem of constructing a circle pattern is transformed into 
the analytic problem of finding the (euclidean or hyperbolic) radii of the circles, 
which have to satisfy some non-linear equations (closure conditions). These non- 
linear closure conditions turn out to be variational: The functionals 5ouc and S'hyp 
(defined in sections 12.31 and \'2 . 4() are functions of the radii, and their critical points 
are the solutions of the closure conditions. The functionals S'euc, 'S'hyp are convex 
(except for scale- invariance in the euclidean case). This implies the uniqueness 
claims of theorem ll.81 fscction 12.511 . The existence claim is more difficult. We have 
to show that the functionals tend to infinity if some radii go to zero or infinity. In 
section 1^71 we show that this is the case if a 'coherent angle system' exists. This 
is a function on the oriented edges, which satisfies a system of linear equations 
and inequalities. The existence problem for circle patterns is thus reduced to the 
feasibility problem of a linear program. In section El^ we prove the existence of a 
coherent angle system if the conditions of theorem ll.8l are satisfied. This is done by 
interpreting a coherent angle system as a feasible fiow in a network (with capacity 
bounds on the branches) and invoking the feasible flow theorem. 

In section 1^21 theorem 1 1.71 for circle patterns in hyperbolic surfaces is deduced 
from theorem 11.81 using methods of combinatorial topology. The basis of this de- 
duction is lemma [2. 141 We present a self-contained proof of it in appendix IdI 

In section I2.1UI the analogous theorem 11.51 for circle patterns in the sphere is 
deduced from theorem ll.71 First, the problem is transferred to the plane by stere- 
ographic projection. Then the proof proceeds in a similar way as in the hyperbolic 
case. 

1.3. Variational principles 

In chapter |21 we define the functionals S'euc, S'hyp, and S'sph for euclidean, 
hyperbolic, and spherical circle patterns. The functional S'sph is not convex. Thus, 
we cannot use it to prove existence and uniqueness theorems. The variables are (up 
to a coordinate transformation) the (euclidean, hyperbolic, or spherical) radii of the 
circles. A Legendrc transformation of these functionals fsection f3.1|l leads to a new 
variational principle involving one new functional S for all geometries (euclidean, 
hyperbolic, and spherical). The variables of S are certain angles; and the variation 
is constrained to coherent angle systems. Depending on whether the constraint 
involves euclidean, hyperbolic, or spherical coherent angle systems (section 12.711 . 
the critical points correspond to euclidean, hyperbolic, or spherical circle patterns. 

Colin de Verdiere first used a variational principle to prove existence and 
uniqueness for circle packings |16j . He constructs two functionals, one for the 
euclidean case, one for the hyperbolic case. The variables are the radii of the cir- 
cles. Critical points correspond to circle packings. Explicit formulas are given only 
for the derivatives of the functionals, not for the functionals themselves. In sec- 
tion E21 

we derive Colin de Verdiere's functionals from our functionals S'euc and 
S'hyp. In particular, this effects the integration of Colin de Verdiere's differential 
formulas. 

Apparently, Bragger |12| had already tried to integrate Colin de Verdiere's for- 
mulas. He derives a new variational principle for euclidean circle packings. The 
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variables of his functional are certain angles, and the variation is constrained to co- 
herent angle systems. This functional turns out to be related to Colin de Verdiere's 
functional by a Legendre transformation. In section 13.41 we derive it from our 
functional S. 

Rivin's functional for euclidean circle patterns |34| is also derived from S {sec- 
tion I3.5|l . It is less general, because the cell decomposition is assumed to be a 
triangulation, and there can be no curvature in the centers of circles. 

Leibon |27| . |28| derived a functional for hyperbolic circle patterns which can 
be seen as a hyperbolic version of Rivin's functional fsection l3.6|l . It is therefore 
natural to expect that Leibon's functional can be derived from S as well. However, 
this is not the case. The Legendre dual of Leibon's functional (section I3.7|l is not 
•Shyp, but a new functional. Unfortunately, we cannot present an explicit formula 
for this functional. In section 14.71 we derive yet another functional, from which 
both S and Leibon's functional can be derived. 

At least since Bragger |12| . there was an awareness of the fact that the circle 
packing functionals have something to do with the volume of hyperbolic 3-simplices. 
Chapter 2] deals with the connection between the circle pattern functionals and the 
volumes of hyperbolic polyhedra. Schlafli's differential volume formula turns out 
to be the unifying principle behind all circle pattern functionals. This geometric 
approach is essential for the construction of the common ancestor of S and Leibon's 
functional. 

Thurston's method to construct circle patterns |45| (implemented in Stephen- 
son's program CirclePack \17\ ) involves iteratively adjusting the radius of each 
circle so that the neighboring circles fit around. This is equivalent to minimizing 
our functionals (S'cuc, •S'hyp) iteratively in each coordinate direction. 

1.4. Open questions 

There is a functional for Thurston type circle patterns (at least in the euclidean 
case), its variables being the radii of the circles. In fact, Chow and Luo |14) show 
that the corresponding closure conditions are variational. Can an explicit formula 
be derived? (For Thurston type circle patterns with "holes" , we derive a functional 
in section rOl ) 

One may also consider Thurston type circle patterns with non-intersecting 
circles. Instead of the intersection angle, one prescribes the inversive distance 
(an imaginary intersection angle) between neighboring circles (see Bowers and 
Hurdal |10j ). Is there a functional for such circle patterns, and can one write 
an explicit formula? Can this approach be used to prove existence and unique- 
ness theorems? These questions are interesting, because inversivc-distance circle 
patterns may be the key to 'conformally parametrized' polyhedral surfaces in E,'^. 

Even though the spherical functional is not convex, may it be used to prove 
existence and (M6bius-)uniqueness theorems for branched circle patterns in the 
sphere? (See also Bowers and Stephenson Can the representation theorem |l.3l 

be generalized for star-polyhedra? This question is interesting because branched 
circle patterns in the sphere can be used to construct 'discrete minimal surfaces' 

Rodin and Sullivan showed that circle packings approximate conformal map- 
pings |37| . Schramm proved a similar result for 'circle patterns with the com- 
binatorics of the square grid' \4:1\ . There have been numerous refinements (for 
example, the proof of 'C°°-convergence' by He and Schramm \21\ ). but all conver- 
gence results deal with circle patterns with the topology of the disc and regular 
combinatorics (square grid or hexagonal). Can the variational approach help in 
proving convergence results for circle patterns with non-trivial topology and irreg- 
ular combinatorics? 
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CHAPTER 2 



The functionals. Proof of the existence and 
uniqueness theorems 

2.1. Quad graphs and an alternative definition for Delaunay circle 

patterns 

A 'quad graph' (the term was coined by Bobenko and Suris jS]) is a cell decom- 
position of a surface such that the faces are quadrilaterals. We also demand that the 
vertices are bicolored. On the other hand, we allow identifications on the boundary 
of a face. For example, figure ITTl (left) shows a quad graph decomposition of the 
sphere with only one 'quadrilateral'. To put is more precisely: 

Definition. A quad graph is a cell decomposition of a surface, such that each 
closed face is the image of a quadrilateral under a cellular map which immerses the 
open cells, and the vertices are colored black and white so that each edge is incident 
with one white and one black vertex. 

From any cell decomposition S of a surface one obtains a quad graph Q such 
that the correspondence between elements of E and elements of Q is as follows: 



s 


Q 


vertices 


black vertices 


faces 


white vertices 


interior edges 


quadrilaterals 



Figure 1231 (right) shows an example which should make the construction clear. 
This construction may be reversed, such that from every quad graph one obtains 
a cell decomposition. (The reverse construction is not quite unique in the case of 




Figure 2.1. Left: The smallest quad graph decomposition of the sphere. 
Right: The quad graph corresponding to the cell decomposition of figure 11.51 (dot- 
ted). Note that boundary edges do not correspond to quadrilaterals. 
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Figure 2.2. Not a Delaunay type circle pattern by the definition of section IlTI 

surfaces with boundary, because one is free to insert any number boundary edges. 
Fortunately, boundary edges play no role in our treatment of circle patterns.) 

The following simple definition of Delaunay type circle patterns in terms of 
quad graphs is a bit more general than the one in section Fl. II 

Definition. A (generalized) Delaunay type circle pattern is a quad graph in 
a constant curvature surface, possibly with conc-likc singularities in the vertices, 
such that the edges are geodesic and all edges incident with the same white vertex 
have the same length. 

This definitions allows for configurations as shown in figure 12.21 The corre- 
sponding cell decomposition has a digon corresponding to the white vertex in the 
middle. 

2.2. Analytic formulation of the circle pattern problem; euclidean case 

Consider the following euclidean circle pattern problem: For a given finite cell 
decomposition S of a compact surface with or without boundary, a given angle 9e 
with < 6*6 < TT for each interior edge e, and a given angle for each face /, 
construct a euclidean Delaunay type circle pattern (as defined in section ri.l|l with 
cell decomposition E, intersection angles given by 0, and cone angles and Neumann 
boundary conditions given by $. 

We will reduce this problem to solving a set of nonlinear equations for the radii 
of the circles. The following lemma is the basis for this. 

Lemma 2.1. Let T, be a cell decomposition of a compact surface, possibly with 
boundary. Let 6 : Eint ^ (0, tt) be a function on the set Eint of interior edges, 
and r : F —I- (0, oo) be a function on the set F of faces. Then there exists a unique 
euclidean circle pattern with cone-like singularities in the vertices and in the centers 
of circles such that the corresponding cell decomposition is S, the intersection angles 
are given by 6 and the radii are given by r. 

The cone angle Qy in a vertex v is given by 0„ = "^Oe, where the sum is taken 
over all edges e around v. The cone angle in the center of an interior face fj (or 
the boundary angle for a boundary face ) is 



where the sum is taken over all interior edges e between the face fj and its neigh- 
bors fk . 

Proof. Given the cell decomposition, intersection angles, and radii, one con- 
structs the corresponding circle pattern as follows. For each interior edge e with 
faces fj and fk on cither side, construct a euclidean kite shaped quadrilateral with 




(2.1) 



2.2. ANALYTIC FORMULATION; EUCLIDEAN CASE 
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-log^^-^. (2.2) 



side lengths rf. and r/^ and angle 9^ as in figure ITSI (left). Glue these quadrilat- 
erals together to obtain a flat surface with cone-like singularities, and, in fact, the 
desired circle pattern. The uniqueness claim is obvious, as is the claim about 6„. 
For each oriented edge e, let (pg be half the angle covered by e as seen from the 
center of the circle on the left side of e. Sec figure IT^ (left). Now, 

where the sum is taken over all oriented edges in oriented the boundary of /. and 

llog"^-'"^-^V" 
2« -^A-e*' 

(The argument of a non-zero complex number z is argz — ^ log |.) Equation (|2.1|l 

follows. ■ 

It is convenient to introduce the logarithmic radii 

p^logr (2.3) 
as variables. Then, equation 1)2. 2() may be rewritten as 

^e = feAPh- Pf,), (2.4) 

where, for < 6 < tt, the function : R ^ R is defined by 

1 1 — e^~*^ 
/«(^)^=^l°gT^^' (2-5) 

and the branch of the logarithm is chosen such that 

< feix) < TT. 

In the following lemma, we list a few properties of the function fe{x) for reference. 
Lemma 2.2. 

(i) The remaining angles of a triangle with sides 1 and and with enclosed 



angle 9 are f0{x) and fg{—x), as shown in fiaure XzlA (right), 
(ii) The derivative of fg{x) is 



- 2(coshr-cosh.) > 0' (2-^) 



so feix) is strictly increasing. 
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(in) The Junction fg{x) satisfies the functional equation 

feix) + fei^x) = 7T - 9. (2.7) 
(iv) The limiting values of felx) are 

lim feix) and lim fe{x) = tt - 6*, (2.8) 



.T^ — OO 



(v) For < y < 7T — 6, the inverse function is 



f^\y) ^log-^^ (2.9) 
sm[y + 0) 

(vi) The integral of fg{x) is 

Fe{x):= f fg{0di = lu,U2{e^+^\ (2.10) 

J —oo 

where Li2(z) is the dilogarithm function; see avvendixJBl 

By lemma 12.11 the euclidean circle pattern problem is equivalent to the non- 
linear equations H2.11() below. 



Lemma 2.3. Given a cell decomposition Yi of a compact surface with or without 
boundary, an angle Oe with < 9^ < tt for each interior edge e, and an angle $ / 
for each face f . Suppose r e and p g are related by equation (|2.3(l . Then 
the following statements (i) and (ii) are equivalent: 

(i) There is a euclidean circle pattern with radii rf, intersection angles 9e and 
cone/boundary angles $/. 

(ii) For each face f £ F , 

^f^2j2feAPf.~Pf)=0, (2.11) 
where the sum is taken over all interior edges e between the face f and its neighbors 

fk- 

2.3. The euclidean circle pattern functional 

The euclidean circle pattern functional defined below is a function of the loga- 
rithmic radii pf. Equations (12.11(1 are the conditions for a critical point. 

Definition. The euclidean circle pattern functional is the function 

<S'cuc(p) = 

(lmLi2 {eP'--P'^+''^) + lmLi2 {e'f.-P!^+''') - (tt - 0,)(py. + p^j) 

+ Y.^fPf- (2-12) 

0/ 

The first sum is taken over all interior edges e, and fj and fk are the faces on either 
side of e. (The terms are symmetric in fj and fk^ so it does not matter which face 
is considered as fj and which as f^.) The second sum is taken over all faces /. 

Lemma 2.4. A function p G E^ is a critical point of the euclidean circle pattern 
functional Seuc, 'if cind only if it satisfies equations (|2.11|) . The critical points of 
Seuc are therefore in one-to-one correspondence with the solutions of the euclidean 
circle pattern problem. 



E 



2.4. THE HYPERBOLIC FUNCTIONAL 
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Figure 2.4. 



Proof. Using equations (|2.1()(l and (|2.7(l . one obtains 



dpf . 

where the sum is taken over all edges e between the face / and its neighbors fk ■ B 

2.4. The hyperbolic circle pattern functional 

This case is treated in the same fashion as the euclidean case. Of course, the 
trigonometric relations are different: 

Lemma 2.5. Suppose ri and r2 are two sides of a hyperbolic triangle, the en- 
closed angle between them is 6, and the remaining angles are ipi and (p2, as shown 
in figure \2.4\ Then 

fi ^ fe{p2 - Pi) - fe{p2 + Pi), (2.13) 
where fe{x) is defined by equation H2.5|l . and 

/9 = logtanh^. (2.14) 

The inverse relation is 

P^ = - log = ^ = ^ (2 15) 



1 

log 

sm r ■ ' ^- Jsm^ ^ 

where 6* = n — 6, ipi^2 > 0, and tpi + ip2 < 0* . 

Equation H2.13() is derived in appendix ^ Equation (|2.15() follows from H2.13|l 
(and the corresponding equation for (^2) by a straightforward calculation using 
equations (|2.9() and l|2.7(l . Note that positive radii r correspond to negative p. 

In the hyperbolic case, the new variables p are given by equation H2.14|l . Instead 
of equation H2.4|) . one has 

= feAPh - Pf,) - feAPh + Pf,)- (2-16) 
and the nonlinear equations for the variables pf are 

- 2 E - Pf) - foAPf. + Pf)) = 0, (2.17) 

where the sum is taken over all interior edges e between the face / and its neighbors 

fk- 

Lemma 2.6. If p & IR^ is a solution of the equations H2.17|l . then pf < 
for all f £ F . The solutions of the equations (|2.17|l are therefore in one-to-one 
correspondence with the solutions of the hyperbolic circle pattern problem. 
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Proof. Since, by equation 1)2. 6|l . the function fe{x) is strictly increasing for 

o<e<n, 

foAPfk + Pf) - foAPh - Pf) > 0. if Pf > 0- 
Since $y > by assumption, the left hand side of equation H2.17|l is positive if 
Pf>0. ■ 



Definition. The hyperbolic circle pattern functional is the function 

'hyp '■ 



Shv-n ■ ^ E- 



5i,yp(p) = J2 (l"^Li2 (6"^'^-"^.+"'^) + ImLi2 {e''f.~'^f^+^'') 

+ ImLi2 (6"^.+"^'=+'''=) +ImLi2 {e-Pf^~P'^+''^)^ (2. 18) 

of 

The first sum is taken over all interior edges e, and fj and fk are the faces on either 
side of e. (The terms are symmetric in fj and fk, so it does not matter which face 
is considered as fj and which as fk-) The second sum is taken over all faces /. 

Lemma 2.7. A function p e is a critical point of the hyperbolic circle 
pattern functional Shyp, if and only if p satisfies equations H2.17|l . In that case, p 
is negative. The critical points of Shyp are therefore in one-to-one correspondence 
with the solutions of the hyperbolic circle pattern problem. 



Proof. Similarly as in the euclidean case, one finds that 

- 2 ^ {fg^iPf, - Pf) - feAPk + Pf)) , (2.19) 



dPf . 

fo\ofk 

such that dShyp — 0, if and only if equations (|2.17|1 are satisfied. By lemma ^IM 
p < follows. ■ 

For future reference, we note that equation (|2.19(l and the proof of lemma IT?)! 
imply 



dpf 



> 0, if Pf > 0. (2.20) 



2.5. Convexity of the euclidean and hyperbolic functionals. Proof of 
the uniqueness claims of theorem 11.81 

Lemma 2.8. // a euclidean circle pattern with data S, 9, $ exists, then the 
euclidean functional is scale- invariant; Multiplying all radii r with the same positive 
factor ( equivalently, adding the same constant to all p ) does not change its value. 

Proof. Let If G be the function which is 1 on every face f G F. Equa- 
tion H2.12fl implies 

^euc(p + /llF) = 5'euc(p) + /l(5I*/-2 J2 ('^ " ^e)) , 

feF eGErnt 

where Eint is the set of interior edges. Clearly, the functional can have a critical 
point only if the coefficient of h vanishes. In this case, the functional is scale 
invariant. ■ 



2.6. THE SPHERICAL FUNCTIONAL 
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If the euclidean functional is scale invariant, one may restrict the search for 
critical points to the subspace 

C/ = {peM^|^p/ = 0}. (2.21) 

Lemma 2.9. The euclidean functional Seuc is strictly convex on the subspace 
U . The hyperbolic functional Shyp is strictly convex on R^. 

Proof. By a straightforward calculation, one finds that the second derivative 
of the euclidean functional is the quadratic form 

^ cosh(,,.-"';)-cos^. 

where the sum is taken over all interior edges e, and fj and fk are the faces on 
either side. Since it is (quietly) assumed that the surface is connected, the second 
derivative is positive unless dpfj = dpf,, for all fj,fk G F- Hence it is positive 
definite on U. 

For the hyperbolic functional, one obtains 

^'-^ ^ Cosh(,,.-p';)-cos.. 



— — (dp/^ +dpf^} 



cosh(p/^ +P/J - CO 

which is positive definite on M^. ■ 
This proves the uniqueness claims of theorem 11.81 

2.6. The spherical circle pattern functional 

Like in the euclidean and hyperbolic cases, there is a functional for spherical 
circle patterns whose critical points correspond to solutions of the circle pattern 
problem. 

Lemma 2.10. Suppose ri and r^ are two sides of a spherical triangle (with 
< ri_2 < T^), the included angle between them is 9, and the remaining angles are 
ipi and if2, CLS shown in figure\2.4\ Then 



VI = fe{p2 - Pi) + f.-B{p2 + Pi), (2.22) 
where fg{x) is defined by equation H2.5|l . and 

p = logtan^. (2.23) 

The inverse relation is 

/oi = - log ^ = ^ = ^ (2 24) 

' 2 S 3in(^: + ^^)sin(^:^ + ^^) ' ^ ■ ^ 

where 6* ~ n — 6, ipi^2 > 0, and 6* < + ip2 < in — 9* . 

Equation H2.22(l is derived in appendix El Equation (|2.24|) follows from p. 22(1 
(and the corresponding equation for (p2) by a straightforward calculation using 
equations (|2.9|) and 12. 7|) . 

To construct the spherical circle pattern functional, one proceeds like in the 
euclidean and hyperbolic cases fsections l2.2H2.4|l . In this case, the new variables pf 
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are given by equation H2.23() . There is a one-to-one correspondence between radii r 
with < r < TT, and p G M. Instead of equations (|2.4() or ()2.16(l . one has 

= /e. {Pf, ~ Pf, ) + fet {Pf, + P/, ), (2.25) 
where 0* — it — 9. Consequently, the nonhnear equations for the variables pf are 

- 2 ^ {feSPf, - Pf) + fotiPf, + Pf)) = 0, (2.26) 

where the sum is taken over all interior edges e around /, and fk is the face on the 
other side of e. 

Definition. The spherical circle pattern functional is the function 

^sph '■ E-^ — > E- 

SsM = E (lniLi2 {eP^k-Pf.+^'^) + ImLi2 (e"/. "''/^ +'«^) 

- ImLia (e''/.+''/.+*(— ''^)) - ImLia (e,-''^. -''^'^+*("-^^)) 

- -^{Pf, + PfS) 
+ Y.^fPf- 

of 

(2.27) 

The first sum is taken over all interior edges e, and fj and fk are the faces on either 
side of e. (The terms are symmetric in fj and /fe, so it does not matter which face 
is considered as fj and which as fk-) The second sum is taken over all faces /. 

Lemma 2.11. A function p G is a critical point of the spherical circle pattern 
functional Ssph, */ if^d only if p satisfies equations (|2.26|l . The critical points of 
Ssph are therefore in one-to-one correspondence with the solutions of the spherical 
circle pattern problem. 

This proposition is proved in the same way as the lemmas 12.41 and 12.71 
The spherical circle pattern functional is not convex. By a straightforward 
calculation, one finds that 

^-'^= ^ (cosh(p,. -"';)- COS.. (^^^^-^^^^)' 
fj°l°fk 

~ cosHpJt)^coseJ '^^^^'^f^)')- ^'-''^ 

This quadratic form is negative for the tangent vector If E , which has a 1 in 
every component. Hence, the negative index is at least one. 

Remark. This has a geometric explanation. Consider a circle pattern in the 
sphere. Focus on one flower: a central circle with its neighbors. The neighbors 
nicely fit around the central circle. Now decrease the radii of all circles by the same 
factor. The effect is the same as increasing the radius of the sphere by that factor. 
This makes the sphere flatter. The neighbors will not fit around the central circle 
anymore, but there will be a gap. To adjust the radius of the central circle so that 
the neighbors fit around, one would make it even smaller. 
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Nevertheless, in numerical experiments, the spherical functional has been used 
with amazing success to construct circle patterns in the sphere. This following 
method was used. Consider the reduced functional 

5,ph(p) - max Ssphip + t If). (2.29) 

Clearly, S'sph(p) is invariant under a shift p p + tip. To solve a spherical cir- 
cle pattern problem, minimize 5'sph(p) on the set U defined by equation H2.21|l . 
Equivalently, minimize Sgphip) under the constraint pf = for some fixed f Cz F. 

The numerical evidence suggests that this method works whenever the circle 
pattern in question exists. In particular, this method can be used to construct 
branched circle patterns in the sphere. In general, such patterns cannot be con- 
structed using the euclidean functional after a stereographic projection. 

Wc return to the definition of the reduced functional Ssphip), equation H2.29|l . 
The following proposition provides a geometric interpretation for the condition 

d 



^^dSsphip + tli 



= 0. 

=0 



Consider the spherical circle pattern problem with data E, 6* S (0,7r)-^, and <!> S 
(0, oo)^. For the sake of simplicity, assume that S is a decomposition of a surface 
without boundary. Suppose there exists a corresponding circle pattern in a surface 
M with constant curvature 1 and cone-like singularities. Then, due to the Gauss- 
Bonnet theorem, the total area A of M is determined by the formula: 

A + K, + Kf = 2T:i\F\ - \E\ + \V\), (2.30) 

where Ky and Kf are as in equations (|1.2|) and and \E\, and \ V\ are the 

numbers of faces, edges and vertices of S, respectively. 

Proposition. Let p G E,^ he arbitrary. Given a cell decomposition Y, of a 
closed compact surface, the intersection angles 9, and radii r ~ 2arctane'', con- 
struct the surface M'^p^ with constant curvature 1 and cone-like singularities by 
gluing together spherical kite-shaped quadrilaterals, as was described (for the eu- 
clidean case) in the proof of lemma \Ej\ Let A^^^ be the total area of AI^p\ Then 

^dSsphip + tlp) ^A-Af") (2.31) 
dt t=o 

where A is determined by equation 1)2.30(1 . 

Proof. Let Ai''^ be the area of the quadrilateral corresponding to an edge e 
as in figure (left). Then 

Ai^^ ^2tp + 2tp' + 29^ - 2-K 

12.251 

= VeXPk - Pj) + VeXPi - Pfc) + 4/(._e,)(Pj -f Pfe) + 20^ - 27r 



^ 4/(^-8e)(Pi +Pfc)■ 



Also, 



and 



^ 29, = 2^\V\ -^K, 
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Finally, 

^ dSsphip + *1f) = ( X! J~) '5'sph(p) 



( - (P, + Pfc) + Vi.-e^) hPj - Pk) - 27r) + ^ 



E (-4/(.-e.)(p,+Pfc)+20e-27r)+^$/ 

Equation l|2.30() follows. ■ 

2.7. Coherent angle systems. The existence of circle patterns 

In section l^^Kl the uniqueness of a circle pattern was deduced from the convexity 
of the euclidean and hyperbolic functionals. This section and the next one are 
devoted to the existence part of theorem 11.81 To establish that the euclidean 
functional attains a minimum, we will show that 

S'euc(p) — * oo if \\p\\ — > oo in U, 

where U £ is the subspace defined by equation H2.21|l . To establish that the 
hyperbolic functional attains a minimum, we will to show that 

Shy-pip) — * if IIpII — ^ °o with p < 0. 

Because of the inequality H2.2()|l . this suffices. 

To estimate the functionals from below, one has to compare the sum over 
interior edges with the sum over faces in equations (|2.12|) and (|2.18|) . This is 
achieved with the help of a so called 'coherent angle system'. In this section, we 
prove that the functionals have minima if and only if coherent angle systems exist. 
In section 12.81 we will show that the conditions of theorem 11.81 are necessary and 
sufficient for the existence of a coherent angle system. 

Coherent angle systems also play an important role in chapter |3 when we 
derive other variational principles by Legcndre transformations. Spherical coherent 
angle systems are defined below, but not used until chapter |3| 

Let E be the set of oriented edges. For an oriented edge e £ E, denote by 
— e £ E the edge with the opposite orientation, and by e the corresponding non- 
oriented edge. 

Definition. A euclidean coherent angle system is a function (p £ M^'"' on the 
set Eint of interior oriented edges which satisfies the following two conditions. 

(i) For all oriented edges e £ Eint, 

ipg > and ipg + ip^g = TT — 6'e. 

(ii) For all faces f £ F, 

where the sum is taken over all oriented interior edges e in the oriented boundary 
of/. 

A hyperbolic coherent angle system satisfies 
(i' ) For all oriented edges e G E, 

ipg > and ipg + p-g < IT — 9e- 

and condition (ii) above. 

A spherical coherent angle system satisfies 
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(i" ) For all oriented edges eeE, 

< < n, n — 6e < fs + f-s < + ^e, and ^^pg — ip-s\ < n — 6e, 

and condition (ii) above. 

(Note that the exterior angles of a spherical triangle satisfy the triangle in- 
equalities.) The following lemma reduces the question of existence of a (euclidean 
or hyperbolic) circle pattern to the question of existence of a coherent angle system. 

Lemma 2.12. The functional Seuc (Shyp) has a critical point, if and only if a 
euclidean (hyperbolic) coherent angle system exists. 

Proof. If the functional Scuc (•S'hyp) has a critical point p, then equation (|2.4(l 
(equation H2.16|l ) yields a coherent angle system. It is left to show that, conversely, 
the existence of a coherent angle system implies the existence of a critical point. 

Consider the euclidean case. Suppose a euclidean coherent angle system (p 
exists. This implies 

$; = 2 ^ (it -9,). 



/GF eeF,„ 



Hence, the functional S'cuc is scale invariant. (See the proof of lemma l2.8l ) We will 
show that S'cuc (p) — > oo if p ^ oo in the subspace U defined in eauation l2.21l More 
precisely, we will show that for p £ U, 

Scuc{p) > 2 min ipg max \pf\. (2.32) 

The functional S'cuc must therefore attain a minimum, which is a critical point. 
For X e M and < 6 < n, 

ImLi2(e^+^'') + ImLi2(e-^+'^) > {tt - 9) \x\, 

and hence, by equation H2.12|l . 

Seuc(p) > -2 ^ (tt - 9e) min (p/, , . ) + ^ $/p/, 

eeE feF 

where the sum is taken over the unoriented interior edges e, and fj and fk are the 
faces on either side of e. Now, we use the coherent angle system ip to merge the 
two sums. Because 

feF eeE,„t 

one obtains 

Scuc(p) > 2 ^ min {pg, p^g) \pf, - P/, |- 

eeE,„, 

Since we assume the cellular surface to be connected, we get 

S'cuc(p) > 2 min pg( max pf — minpf), 
eeE.„, f&F feF'^" 

and from this the estimate H2.32|l . 

The hyperbolic case is similar. One shows that, if all p/ < 0, 

S'hyp(p) > 2 min \pg + p^g ~ {-k — 9e)\ maxlp/l. 

eSF.nt feF 
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Figure 2.5. The network {N,X). Only a few of the branches and capacity inter- 
vals are shown. 

2.8. Conclusion of the proof of theorem 11.81 

With this section we complete the proof of thcorcm ll.8l In scction [T5l we have 
shown the uniqueness claim. By lemma 121 the circle patterns exist, if and only 
if a coherent angle system exists. All that is left to show is the following lemma. 

Lemma 2.13. A euclidean/hyperbolic coherent angle system exists if and only 
if the conditions of theorem \l.Si hold. 

The rest of this section is devoted to the proof of lemma [^.131 It is easy to see 
that these conditions are necessary. To prove that they are sufficient, we apply the 
feasible flow theorem of network theory. Let {N, X) be a network (i.e. a directed 
graph), where A'' is the set of nodes and X is the set of branches. For any subset 
N' C N let ex{N') be the set of branches having their initial node in N' but not 
their terminal node. Let in{N') be the set of branches having their terminal node in 
A'^' but not their initial node. Assume that there is a lower capacity bound and an 
upper capacity bound associated with each branch x, with — oo < < < oo. 

Definition. A feasible flow is a function (p e R-^ , such that Kirchoff's current 
law is satisfied, that is, for each n G N, 

xGex{{n}) x^in{{n}) 

and fla: < </3x < bx for all branches x €z X. 

Feasible Flow Theorem. A feasible flow exists if and only if for every 
nonempty subset N' C N of nodes with N' ^ N, 

'Y ^ flj;. 

x£ex{N') x£in{N') 

A proof is given by Ford and Fulkerson |18[ ch. II, §3]. (Ford and Fulkerson 
assume the capacity bounds to be non-negative, but this is not essential.) 

To prove lemma E. 131 in the euclidean case, consider the following network; see 
figure The nodes are all faces and non-oriented interior edges of the cellular 
surface, and one further node that we denote by H: N = F U E U {K}. There is 
a branch in X going from H to each face f G F with capacity interval [^$/,oo). 
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From each face / there is a branch in X going to the non-oriented interior edges 
of the boundary of / with capacity interval [e, oo), where e > will be determined 
later. Finally there is a branch in X going from each non-oriented edge e G i? to 
Kl with capacity (— cxd, tt — 0e]. 

Assume the conditions of theorem ll.8l are fulfilled. A feasible flow in the network 
yields a coherent angle system. Indeed, since we have equality in H1.5|l if F' = F ^ 
Kirchoff's current law at M implies that the flow into each face / is and the 
flow out of each edge e \s tt — 9e- It follows that the flow in the branches from F to 
E constitutes a coherent angle system. 

We need to show that the condition of the feasible flow theorem is satisfied. 
Suppose N' is a nonempty proper subset of N. Let F' = N'dF and E' ~ N'dEint- 

Consider first the case that Kl e N\ which is the easy one. Since A^' is a proper 
subset of A^ there is a face f £ F or an edge e £ E which is not in A^'. In the first 
case there is a branch out of A^' with infinite upper capacity bound. In the second 
case there is a branch into A^' with negative infinite lower capacity bound. Either 
way, the condition of the feasible flow theorem is trivially fulfilled. 

Now consider the case that M ^ N' . We may assume that for each face f (1 F' , 
the interior edges in the boundary of E are contained in E' . Otherwise, there would 
be branches out of A^' with infinite upper capacity bound. For subsets A,BcN 
denote hy A ^ B the set of branches in X having initial node in A and terminal 
node in B. Then the condition of the feasible fiow theorem is equivalent to 

^ l<^j + e\F\F'^E'\< 



feF' eeE' 



It is fulfilled if we choose 



1 

e < — ; — r mm 
2E F 



'n( E i^-ee)-J2l'^X (2.33) 

\eeE'{F') feF' J 

where the minimum is taken over all proper nonempty subsets F' of F and E'{F') 
is the set of all non-oriented edges incident with a face in F' . The minimum is 
greater than zero because strict inequality holds in (|1.5|l if F' is a proper subset of 
F. 

In the hyperbolic case, the proof is only a little bit more complicated. In the 
network, the flow in the branches going from B to a face / must constrained to be 
exactly 5$/; and the capacity interval of branches going from an edge e to K has 
to be changed to (—00, 9^ — e]. 

2.9. Proof of theorem 11.71 

The necessity of the condition of theorem 11.71 can be seen from the following 
geometrical argument. Figure 12.61 shows a closed path in the Poincarc dual cell 
decomposition S*, which cuts out a disc. By the Gauss-Bonnet theorem, 

E(7r - aj - Pj - 7j) = 27r - eA, 

j 

where e e {0, —1} is the curvature of the surface and A is the area enclosed by the 
path. But 

TT - a^- - I3j = 6j - eAj, 
where Aj is the are of triangle j. Hence, 



3 
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Figure 2.6. The circles and quadrilaterals corresponding to a closed path in the 
Poincare dual (dashed), which cuts out a disc. 

Inequality H1.4|l follows, with equality holding if and only if all are zero and 
e{A = 0- This is the case if and only if the path contains only a single 

vertex of E in its interior. 

We will now prove the sufficiency of the condition. The principle tool in the 
proof is the following lemma. 

Lemma 2.14. Suppose a nonempty 1- dimensional cell complex (a graph) T = 
{Vi,Ei) is cellularly embedded in a cell decomposition E = [V,E,F) of a closed 
compact surface. Suppose T separates S into r contiguous regions. Then 

r 

r-\E^\ + \V^\ = \F\ - \E\ + \V\+Y, h„ (2.35) 

where hj is the dimension of the first Wi2-homology group of the j*^ region. 

If all regions are simply connected, then all hj = 0, and equation H2.35(l states 
that the Euler characteristic is invariant under subdivisions. A multiply connected 
region can be turned into a simply connected one by adding hj cuts to F which do 
not separate the region. Each cut reduces the left hand side of equation (|2.35|) by 
one. For a more formal proof of lemma 12. 141 see Giblin |20l eh. 9] or appendix IdI 

Suppose the condition of theorem 1 1 . 71 holds . We will deduce the condition of 
theorem 1 1.81 Thus, let F' C F be a nonempty subset of the set F of faces, and let 
E' C E he the set of edges, which are incident with any face in F' . We have to 
show inequality (|1.5|) where all ~ 2tt. That is, we have to show 

291 > 27r|F'|, (2.36) 

eeE' 

where equality holds if and only if the genus g = 1 and F' = F. 
Let 

F" ^F\ F' and E" ^ E\ E' 
be the complements of F' and E' . Identify the faces, edges, and vertices of the cell 
decomposition S with the corresponding vertices, edges, and faces of its Poincare 
dual S*. Consider the 1-dimensional subcomplex F = {F",E") of S* which con- 
sists of the vertex set F" and edge set E". (It is easy to see that it is indeed a 
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subcomplex.) First, suppose that F' ^ F. Then the graph F is not empty. By 
lemma I^HI apphed to F in S*, 



\E"\ + \F"\ = \V\^\E\ + \F\+Y.h, 



or, equivalently. 



Now, 



\F'\^\E'\ + \V\^Y.^l~h,). (2.37) 



^ 2ei=2n\E'\-Y^ 29,. 

2'k\E'\ - ^ 26ie + 26, 



Because, by the condition of theorem 1 1.71 sums to 2tt aroimd each vertex, 

eG-E 

and hence 

^ 291 = 2^{\E'\ - \V\) + 29,. (2.38) 
With equation H2.37|l . we get 

r 

Y 2^ - 2^l^'l = E - 27r^(l - /i,). (2.39) 

eG£" eeB" j=l 

Let e{, . . . ,el. be the boundary path of the j"* region into which E* decomposes if 
cut along the edges in E" . (The same edge may appear twice in one boundary path.) 
Each edge of E" appears exactly twice in any of the boundary paths. Therefore, 

r bj 
ee_E" j = l i=l 

and, with equation H2.39|) . 

r bj 

J2 29: - 27r\F'\ = E (E^e^ - 2^(1 " ^"^^ 

ee-E' i=l i=l 

In any case, 

J20,.^27:(l~h,)>0. (2.40) 

For if the j*'' region is not a disc, then hj > 1, and (|2.40() follows because 9 > 0. 
If the J*'' region is a disc, then hj = and (|2.4U|) follows by the condition of 
theorem 1 1.71 We have shown that inequality H2.36|l holds strictly, if F' ^ F. 
Now assume F' = F, and hence E' = E. Like equation (|2.38(l . we get 

Y20:=27t{\E\-\V\). 

eGE 

Therefore, since 2 - 2g = |F| - \E\ + \V\, 

Y20:^27r{\F\~2 + 2g). 
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Figure 2.7. Left: The regular cubic pattern after stereographic projection to the 
plane. Right: The Neumann boundary condition <I>j for the new boundary circles. 

Thus, inequality H2.36|l holds strictly, except if = |F| and g = 1. This completes 
the proof of theorem 11.71 



A circle pattern in the sphere may be projected stereographically to the plane, 
choosing some vertex, Voo, as the center of projection. (For example, figure \T7\ 
(left) shows the circle pattern combinatorially equivalent to the cube and with 
intersection angles 7r/3 after stereographic projection.) One obtains a circle pattern 
in the plane, in which some circles (those corresponding to faces incident with Voo) 
have degenerated to straight lines. Since stereographic projection is conformal, the 
intersection angles are the same as in the spherical pattern. Furthermore, Mobius- 
equivalent circle patterns in the sphere correspond to patterns in the plane which 
differ by a similarity transformation; that is, provided the same vertex is chosen as 
the center of projection. 

To prove that the condition of theorem 11.51 is necessary for a circle pattern 
to exist, project the pattern to the plane as described above and proceed as in 
section EHl Some circles may have degenerated to straight lines, but equation 1)2.34(1 
holds nonetheless, with e = 0. All jj are zero if the dual path has one finite vertex 
in its interior, or if all circles have degenerated to straight lines. In that case the 
dual path on the sphere encircles the vertex which is the center of projection. 

It is left to show that the condition of theorem 1 1.51 is sufficient. So assume the 
condition holds. The idea is to show that corresponding planar pattern exists using 
theorem 1 1.81 and then project it stereographically to the sphere. 

To show the existence of the planar pattern, first choose a vertex Voo of the cell 
decomposition E of the sphere. Let F be the set of faces of E, and let Foo C F he 
the set of faces which arc incident with v^o ■ Then remove from E the vertex Voo , 
all the faces in Foo, and the edges between them, to obtain a cell complex Eq with 
face set Fq = F \ Foo ■ Because E is a strongly regular cell decomposition of the 
sphere, Sq is a cell decomposition of the closed disc. Hence, theorem 11.81 mav be 
applied to prove the existence and uniqueness of the planar pattern. (Except in the 
trivial case when Eq has only one face.) The Neumann boundary conditions need 
to be specified. For a boundary face / of Eg, set 



where the sum is taken over all boundary edges e of Eg which arc incident with /. 
See figure (right). For all interior faces /, set 
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(2.41) 



$/ = 27r. 



(2.42) 
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If the conditions of theorem ll.Sl are satisfied, one may construct the corresponding 
planar pattern, add the hnes corresponding to the removed faces, and project to 
the sphere. Hence, we need to show the statements (i) and (ii): 

(i) For a boundary face / of Eqi ^/ as defined in equation (|2.41|) is positive. 

(ii) If Fq e Fq is a nonempty set of faces of Eq, and E'^ is the set of of all 
interior edges of Eq which are incident with any face in fg, then 

^ < ^ , (2.43) 

where equality holds if and only if i^Q = Fq. 

Lemma 2.15. The statements (i) and (ii) above follow from the statement (Hi) 
below. 

( Hi ) If F' is a nonempty subset of \ Fo^ , and E' is the set of all edges of E 
which arc incident with any face in F' , then 

2^1^^'! < J2 2^*(e)' (2-44) 

where equality holds if and only if i^' = F \ Foo ■ 

Proof. Wc deduce statement (i). First let F' = F \ Foo to obtain 

2n\F'\ = 2C- (2.45) 

edges 
of So 

Then let F' = F \ (Foo U {/}), where / is a boimdary face of Eq. Subtract the 
corresponding strict inequality from equation H2.45|l to obtain 

2^>^2C, 

where the sum is taken over all boundary edges of Eq which arc incident with /. 
Hence, assertion (i) is true. 

We deduce statement (ii). If Fq = F' then 

Fq = F' \ {boundary edges of Eq}. 

By the definition of $/, inequality H2.44|l is equivalent to inequality H2.43|l . ■ 

It is left to prove assertion (Hi) under the assumption of the condition of the- 
orem ^3] We proceed in a similar way as in section 

Suppose that F' is a subset of F \ Foo , and F' is the set of all edges of E which 
are incident with any face in F'. Consider F" = F \ F', the complement of F', and 
E" = E\E' , the complement of F'. Consider the Poincare-dual cell decomposition 
E* and, in it, the 1-dimensional subcomplex (or graph) P = (F",F") with vertex 
set F" and edge set F". As for any graph, we have 

|F"|-|^"|-c-n, 

where n is the number of connected components of P and c is the dimension of the 
cycle space. Since the graph is embedded in E*, a cellular decomposition of the 
sphere, we have 

c = r — 1, 

where r is the number of regions into which P separates E* . Since E" contains the 
edges of E incident with Woo, or, dually, the edges of E* in the boundary of Woo, the 
number of regions is at least two. Hence, the boundary of each region is nonempty. 
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By the condition of thcorem ll.5l the sum of 6 over the boundary of each region is 
at least 27r. Sum over all regions to obtain 

27rr < ^ 29 f,. 

eeE" 

Indeed, each edge in E" appears in or 2 boundaries. Equality holds if and only if 
every edge of E" is in the boundary of a region and each boundary is the boundary 
of a single face of E*. Thus, equality holds if and only if E" ^ E or E" is the 
boundary of a single face of E*. This is the case, if and only if F' = (this is ruled 
out by assumption) or F' = F \ F^o ■ 
Thus, we have shown that 

27ri\E"\ - |F"|) < ^ 20e - 2n{n + 1), (2.46) 

eeE" 

with equality if and only if F' = F \ Foo ■ 

Equation H2.38|l from section 12.91 also holds here. Thus, inequality (|2.46(l is 
equivalent to 

-2^(|F| - |F| + |y|) + < ^ 29; 2n{n + 1), 
and, using Euler's formula, 

\F\-\E\^\V\^2, 

equivalent to 

27r|F'| < ^ 261* - 27r(7i - 1). 

eG-E' 

Since n > 1, and n = 1 if F' = F \ Foo, we have deduced the assertion (in) above. 
This completes the proof of theorem ll.5l 



2.11. Proof of theorem ITTsl 

There is a one-to-one correspondence between Delaunay type circle patterns 
in the sphere and polyhedra in hyperbolic 3-space with vertices in the infinite 
boundary. In the Poincare ball model, hyperbolic space corresponds to the interior 
of the unit ball. The unit sphere corresponds to its infinite boundary. Hyperbolic 
planes are represented by spheres that intersect the unit sphere orthogonally. Hence, 
there is a correspondence between circles in the unit sphere and hyperbolic planes. 
Furthermore, the intersection angle of two circles equals the dihedral angle of the 
corresponding planes. The isometrics of hyperbolic space correspond to the Mobius 
transformations of the sphere at infinity. 

2.12. Proof of theorem [TT^ 

Theorem II. 51 has the following corollary. 

Corollary. Let T, be a strongly regular cell decomposition of the sphere, and 
suppose every vertex has n edges. (Because T, is strongly regular, n > 3.) In other 
words, in the Poincare dual decomposition T,* , every boundary of a face has n edges. 
Suppose that every simple closed path in E* which is not the boundary of single face 
is more than n edges long. Then there exists, uniquely up to Mobius transforma- 
tions, a corresponding circle pattern in the sphere with exterior intersection angles 
2TT/n. 
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Figure 2.8. A cellular decomposition (left) and its medial decomposition (right). 

The case n ~ A implies theorem 11.21 Indeed, suppose S is a strongly regular 
cell decomposition of the sphere. In theorem 11.21 circles correspond to faces and 
vertices. To apply the corollary, consider the medial cell decomposition of 
E. The faces of !],„ correspond to the faces and vertices of S and the vertices of 
Em correspond to edges of E. Figure 12.81 shows part of a cell decomposition of 
the sphere (left) and its medial decomposition (right). The dotted lines in the left 
figure represent the edges of the medial decomposition. On the right, the faces of the 
medial decomposition which correspond to vertices in the original decomposition 
are shaded. 

The vertices of the medial decomposition E^ are 4-valent. The assumption 
that E is strongly regular implies that E™ is also strongly regular. It also implies 
that every simple closed path in the Poincare dual E^ which is not the boundary of 
single face is more than 4 edges long. See the remark on pagcQ point (vi). Hence, 
the corollary implies theorem II. 21 

2.13. Proof of theorem [Ql 

The main part of theorem 11.31 follows directly from theorem 11.21 Given the 
cell decomposition E, by theorem 1 1.21 there is a Mobius- unique circle pattern with 
orthogonally intersecting circles corresponding to the faces and vertices of E. The 
circles which correspond to the faces lie in planes which form the sides of the 
polyhedron in question. Indeed, the circles corresponding to neighboring faces of 
E touch, and hence the corresponding planes intersect in a line which touches the 
sphere. Consider the faces of E which are incident with one vertex v of E. The 
corresponding planes go through one point, namely, the apex of the cone touching 
the sphere in the circle corresponding to the vertex v. Thus, one obtains a convex 
polyhedron — provided, that is, that all vertex-circles are smaller than a great circle. 
But one can always achieve this by a suitable Mobius transformation. We will show 
below that, by applying a suitable Mobius transformation, one can always get the 
center of gravity of the intersection points of the orthogonal pattern into the center 
of the sphere. In that position, all circles must be smaller than a great circle, 
because otherwise all intersection points would lie in one hemisphere. 

The construction above is reversible. Given a polyhedron with edges tangent 
to the sphere, one obtains a circle pattern as in theorem 1 1.21 

The dual polyhedron is obtained by interchanging the role of face-circles and 
vertex-circles. 

Any Mobius transformation of the sphere is the restriction of a projective trans- 
formation of the ambient space which maps the sphere onto itself. Conversely, a 
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Figure 2.9. The 'distance' to an infinite point v is measured by cutting off at 
some horosphere through v. Left: Poincare ball model. Right: half-space model. 



projective transformation of the ambient space which maps the unit sphere onto 
itself induces a Mobius transformation of the sphere. Thus, the uniqueness claim 
of theorem II .31 follows form the uniqueness claim of theorem ll.2l 

It is left to show that, by applying a suitable Mobius transformation to the 
orthogonal circle pattern, one can get the center of gravity of the intersection points 
into the center of the sphere. This follows from the following lemma. The remainder 
of this section is devoted to proving it. 

Lemma 2.16. Let vi, . . . ,v„ be n>3 distinct points in the d- dimensional unit 
sphere S'^ C B,''^^. There exits a Mobius transformation T of S'^ , such that 

n 

i=i 

// T is another such Mobius transformation, then T = RT , where R is an isometry 
ofS^. 

The proof relies on the close connection between the Mobius geometry of S'^ 
and the geometry of {d + l)-dimensional hyperbolic space H'^^^ . In the Poincare 
ball model of hyperbolic space, is identified with the unit ball in 11"^+^, its 

infinite boundary is 5''. The isometrics of H'^^'^ extend to Mobius transformations 
of S'^. Conversely, every Mobius transformation of S^^ is the extension of a unique 
isometry of H'^^^ . 

Given n > 3 points wi, . . . , w„ € S"**, we are going show that there is a unique 
point X S H'^^^ such that the sum of the 'distances' to wi, . . . , u„ is minimal. Of 
course, the distance to an infinite point is infinite. The quantity to use is the 
distance to a horosphere through the infinite point. Sec figure 

Definition. For a horosphere h in H''^^, define 

{ — dist(a;, h) if x is inside h, 
a X e h, 

dist(x, h) if X is outside h, 

where dist(.T, h) is the distance from the point x to the horosphere h. 

Suppose V is the infinite point of the horosphere h. Then the shortest path 
from X to h lies on the geodesic connecting x and v. If h' is another horosphere 
through V, then Sh — Shi is constant. If 51 : R ^ is an arc-length parametrized 

geodesic, then Sh o g is a, strictly convex function if v is not an infinite point of the 
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geodesic g. Otherwise, Sh o g(s) = ±(s — sq). These claims are straightforward 
to prove using the Poincare half-space model where v is the infinite point of the 
boundary plane. Also, one finds that 



where hj are horospheres through different infinite points and n > 3. Thus, the 
following definition is proper. 

Definition (and Lemma). Let wi, ...,«„ be n points in the infinite boundary 
of H"^^^, where n > 3. Choose horospheres hi, . . . ,hn through vi, . . . , u„, respec- 
tively. There is a unique point x G H'^'^^ for which X]j=i (x) is minimal. This 
point X does not depend on the choice of horospheres. It is the point of minimal 
distance sum from the infinite points Wi, . . . , u„. 

There seems to be no simpler characterization for the point of minimal distance 
sum. However, it is easy to check whether it is the origin in the Poincare ball model. 

Lemma 2.17. Let vi, ...,?;„ be n > 3 different points in the infinite boundary 
of H'^^^ . In the Poincare ball model, Vj £ S'^ C M.'^'^^ . The origin is the point of 
minimal distance sum, if and only if^ Vj ~ 0. 

Proof. If hj is a horosphere through Vj, then the gradient of at the origin 
is the unit vector (The metric is ds^ = (yt^-^)^ X^'^^j-) ' 

Lemma 12.161 is now almost immediate. Let x be the point of minimal distance 
sum from the wi ,...,«„ in the Poincare ball model. There is a hyperbolic isometry 
T which moves x into the origin. If T is another hyperbolic isometry which moves 
X into the origin, then T = RT, with R is an orthogonal transformation of R'^"'"^. 
Lemma 12.161 follows. 

This concludes the proof of theorem 1 1.31 
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Other variational principles 

3.1. Legendre transformations 

In this section, we derive different variational principles for circle patterns by 
Legendre transformations of the functionals S'ouc, S^ypi and S'sph- We obtain the 
functional S, defined below, which depends not on the (transformed) radii p £ 
but on the angles (p E M,^'"'. According to whether the variation is constrained to 
the space of euclidean, hyperbolic, or spherical coherent angle systems, the critical 
points correspond to circle patterns of the respective geometry. On the euclidean 
and hyperbolic coherent angle systems, the functional S is strictly convex upwards, 
so that there can be only one critical point, which is a maximum. (A function / on 
a convex domain D is called strictly convex upwards ( downwards ), if 

/((I - t).Ti + 1X2) (>) (1 - t)f{xi) + tf{x2) 

for all xi,X2 & D and < t < 1.) The function S{if) is defined in terms of Clausen's 
integral Cl2(a;), see appendix IbI 

Theorem 3.1. Let Eint he the set of oriented interior edges and define 

S : E,^-"' E, 

si.^) = ( {e: + </^e - ^-e) + CI2 [ei - + V-s) 

e: eji — e 

+ ci2(c + + f-s) + ci2(e: - - f-s) - 2 Gh{2e:S) , 

(3.1) 

where the sum is taken over all non-oriented interior edges e, and e, — e are the 
two oppositely oriented edges corresponding to it. 

(i) The function S is strictly convex upwards on the set of all euclidean coherent 
angle systems. 

If the function Seuc attains its minimum at p £ M.^ , then the restriction of S 
to the space of euclidean coherent angle systems attains its maximum at if £ R^"' 
defined by equations H2.4() . 

Conversely, suppose the restriction of S to the space of euclidean coherent angle 
systems attains its maximum at tp £ ]R^. Then the equations 

e 

are compatible. (There is one equation for each oriented interior edge e and fj, 
fk are the faces to its left and right, respectively.) Therefore, they define p £ 
uniquely up to an additive constant. The function Seuc attains its minimum at 
this p. 

If they exist, the extremal values are equal: 

min Seucif) ~ max S{tp) 

pGR,^ ipeuc. coherent 
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(ii) The function S is strictly convex upwards on the set of hyperbolic coherent 
angle systems. 

If the function Shyp attains its minimum at p E E,^, then the restriction of S to 
the space of hyperbolic coherent angle systems attains its maximum at the ip E E,^'"' 
defined by equation (|2.16|l . 

Conversely, suppose the restriction of S to the space of hyperbolic coherent angle 
systems attains its maximum at ip G ]R^'"* . Then the equations 

sin V^f -V-e ) ( 0* ^+^-e \ 

Pf = - log = (3.3) 

' 2 (^:±^|±^) (^:±^|^) 

are compatible and define, therefore, a unique p G R^. (There is one equation for 
each oriented interior edge e, and f is the face on its left side.) The function Shyp 
attains its minimum at this p. 

If they exist, the extremal values are equal: 

min Shyp{p) = max S{lp) 

pelR--^ if^E, hyp. 

coherent 

(Hi) If p G is a critical point of the function Ssph, then (p G E,^'"' defined by 
equations (|2.25() is a critical point of the restriction of S to the space of spherical 
coherent angle .systems. 

Conversely, suppose the restriction of S to the space of spherical coherent angle 
systems is critical at ip G M.^ . Then the equations 

sin f + sin 'Pf +^-e ) 
Pf=^- log ^ (3.4) 



sm 
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are compatible and define, therefore, a unique p G M.^ . (There is one equation for 
each oriented interior edge e, and f is the face on its left side.) The function Sgph 
is critical at this p. 

The values at corresponding critical points are equal: 

S sph{,P critieal) — S (^(p critical) • 

Remark. 1. If is a euclidean coherent angle system, then equation (|3.1(l 
simphfies to 

S{ip) = J2{Ch{2ip,) + Ch{2ip^,) - Cl2(2C)), (3.5) 

where the sum is taken over aU interior non-oriented edges e, and e, — e are the 
corresponding oriented edges. 

2. Equations (|3.3|l and (|3.4(l may be subsumed under the equation 

Zl^p^) Sin (^ (^^-'Pe^+^-e ' 



Pi ^\ log 



sm 



±^PS±iP-e\ ■ (0 +(ps-iP-S 



sm I '-^ 1 sm , 2- 



The rest of this section is devoted to the derivation of the functional S{if) and 
the proof of theorem 13. II 

In classical mechanics, the Legendre transformation translates between the La- 
grangian and Hamiltonian descriptions of a mechanical system. (For a thorough 
treatment, see, for example, Arnold's textbook |3|.) The motion of such a system is 
typically described by functions q{t) which are critical for the functional J L{q, q) dt, 
where L{v, q) is the Lagrangian of the system. The Hamiltonian II(p, q) is obtained 
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by a Legendre transformation with respect to v: H{p,q) = pv — L{v,q), where 
p = dL{v, q)/dv. It turns out that the motion of the system is described by func- 
tions p{t), q{i) which are critical for the functional J [pq — H{p, q)) dt. Here, p and 
q vary independently. 

The construction of S proceeds in a similar way. The role of the position 
variables q is played by the variables p. Whereas q depends continuously on the 
time t, p is a function on the finite set of faces. Instead of integrals, the functionals 
are finite sums. For each functional, euclidean, hyperbolic, and spherical, we will 
define a Lagrangian form first. From it, we will obtain a Hamiltonian form by a 
Legendre transformation. A reduction leads to the functional S in every case. 

First, we recapitulate the definition and basic properties of the Legendre trans- 
formation of a smooth function of one variable. Then we deal with the euclidean, 
hyperbolic, and spherical functionals separately. 

3.1.1. The Legendre transformation of a smooth function of one vari- 
able. The Legendre transformation is defined for convex functions. In general, 
they need not be smooth and may depend on many variables. For the following 
sections, it is sufficient to consider only the simplest case. The generalization to 
many variables, which is needed in section Hl?! is straightforward. 

Definition. Suppose F is a smooth real valued function on some interval in 
R, and F" > (or F" < 0). Then 

y = F'{x) (3.6) 
defines a smooth coordinate transformation. The Legendre transform of F is 

F{y)^xy-F{x), (3.7) 

where x is related to y by (|3.6() . 

The Legendre transformation is an involution: If F is the Legendre transform 
of F and, say, F" > 0, then F" > and F is the Legendre transform of F. Indeed, 

F'iy) = -rV + x- -rF\x) X 
ay dy 

and 

^ dr 1 

F"[y) = ^ 



dy F"{x) 

3.1.2. The eucUdean functionaL Equation (|2.12|l . defining the functional 
Scuc, may be rewritten as 

5cuc(p)= Yl (^e,(p/. -p/,)-(7r-0e)p/J+^$/P/. 

Here, the first sum is taken over all interior oriented edges e, and fj and are 
the faces on the left and right of e, and e is the corresponding non-oriented edge. 
The second sum is taken over all faces /. The function Fg{x) is defined in equa- 
tion 1)2. 10|) . Define the Lagrangian form of the functional by 



iL 



5i 

If, for all oriented interior edges e&E, 
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where fj and fk to are the faces to the the left and to the right of e, respectively, 
then 

Stuci'V^P) = S'cuc(p). 

In general, however, we do not assume that v comes from a p in this way, or even 
that vg = —V-g. 

By equations (|2.6(l and 1)2. 10(1 . Fg{x) > 0. The Lcgendrc transform of Fg{x) is 
Feiy) = Cl2(2y) + 1 Cl2(2y + 29)-^ Ch{2e) , 

where 

This follows from equation l|B.7|l of appendix El By equation l|2.8|l . the domain of 
Fe{y) is 

< y < t: - 6. 
Define the 'Hamiltonian' form of the functional by 

Seuci^, P)= i^siPh - Pf,) - FeAvs) " (tt - 0e)p/J + Yl ^fPf^ 

where the (^-domain D is 

Z) = G E.^"" I < V3e- < TT - 6le for all e e E^] ■ (3.8) 

(This is a superset of the space of euclidean coherent angle systems.) Writing it all 
out, one has 

5fuc(^,P) = E (f^^Ph - Ph) + \{Gh{2^e) - Cl2(2^e-+ 20e) + Cl2(20,)) 



of 



If ip and p are related by equations 1)2.4(1 . then 

Scucif^P) = Scuc{p)- 



Now, Fg > 0. This follows from Fg > and the properties of the Legendre trans- 
formation. (It can also be verified by a direct calculation.) Therefore, S^^A(p,p) is 
strictly convex upwards with respect to (p. 

The motivation for this construction of S^^^ is, that the critical points of S'ouc 
correspond to critical points of S^^^ where p and p may vary independently: 

Lemma 3.1. // {p,p) is a critical point of S^^^, then p is a critical point of 
Seuc- Conversely, suppose p is a critical point of Seuc- Define p by equations (|2.4|l . 
Then {p,p) is a critical point of Sj^^^. At a critical point, Seuc{p) = S^^^p^p). 

Proof. For any interior oriented edge eGE, 

"""-(^^P) = P/fc - Pfj ~ F'eXvs), 



dps 

Hence, all partial derivatives with respect to the variables ipg vanish, if and only if </? 
and p are related by pf^^ - pf. = Fg^pg), or, equivalently, by pg = Fg^pf^ -pjj = 
feAPfk - Pf,)- In this case, S^^^p, p) = Soudp)- ■ 
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Note that S^^^{(p,p) depends linearly on p. Collecting the coefficients of each 
Pf, one finds that, if is a euclidean coherent angle system, then S^^^{tp,p) does 
not depend on p at all. In fact, in that case, 

S^^,i^,p) = Si^). 

To see this, note that if ip is a coherent angle system, then tpg + 9^ = tt — ips- 
Hence Ch{2ips + 29^) = - Ch{2ip_g). Also, 012(26*) = - Cl2(26'*). 

The convexity claim of theorem l3.1l^ z) follows from the convexity of S^^^{(p, p) 
with respect to ip. Equations l|3.2|l are the inverse relations to equations (|2.4() . It is 
left to prove that they are compatible if (/s is a critical point of S under variations in 
the space of coherent angle systems. The rest follows from lemma mi Suppose ip is 
a critical point of S under variations in the space of coherent angle systems. Since 
we are only interested in the restriction of S to the space of euclidean coherent 
angle systems, we may use equation (|3.5|l for S. The partial derivatives are then 

dS . 

- — = -21og(2sm(pg-)- 
dipg 

The tangent space to the space of euclidean coherent angle systems is spanned by 
vectors of the form 



d d 



t-^\dips dip^g/' 

where 7 is some oriented closed path in the Poincare dual cell decomposition. (If 
we add to ipg we have to subtract the same amount from ip-g. But the sum of 
tp around each face has to remain constant, so we step from one face to the next, 
adding and subtracting along some closed path of the Poincare dual.) Since ip is a 
critical point under variations in the space of coherent angle systems, 

E, sin (pg 

^°^in(^. + ^.)=0 

for all closed paths 7 in the Poincare dual. Hence, equations l|3.2|l are compatible. 
This completes the proof of part (i) of theorem 13. II 

3.1.3. The hyperbolic functional. The construction in this case is similar 
to the euclidean case, but a bit more involved. This is due to the fact that the 
nonlinear terms in equation (|2.18|) (the definition of 5hyp(/o)) depend not only on 
the differences pk — pj, but also on the sums pk + pj- In fact, equation H2.18|l may 
be written as 

Shypip) = E i^<^' ^Pf" ~ Pf^ ) + "^o- ^Pf^ + Pf" )) + E 

where 

Gg{x) = Fe{x) + Fg{-x) = ImLia (e"+*^) + ImLia (e-"+'^). 
The first sum is taken over all non-oriented edges e, and fj and are the faces on 
either side of e. The function Gg{x) is even, so that it does not matter which face 
is taken to be fj and which fk- The second sum is taken over all faces /. 

We introduce some notation before we define the Lagrangian form of the func- 
tional. The space K^'"* of functions on the oriented interior edges splits into the 
direct sum of antisymmetric and symmetric functions: 

M^"" = Alt(^,„t) e Sym(^„t), 

where 

Ah(^„t) ^{ve R^"" I v^g = -vg] (3.9) 
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and 

Define the function 

Styp : k\t{E.M) X Sym(^,„t) x ^ M, 

The first sum is taken over all non-oriented interior edges e, and e is one of the 
oriented representatives e. The second sum is taken over all faces /. 
If, for all oriented edges eGE, 

= Ph - Pfj and ws = p/, + p/, , 
where fj and fk are the faces to the the left and to the right of e, respectively, then 

'5'hyp(w,W,p) = 5'hyp(p). 

By a straightforward calculation, 

G'g{x) = fe{x) - fei-x) = 2arctan (tan (^) tanh (|)) , 

and 

= = 2/^(2:) >0. 

The Legendre transform of Gg{x) is 

Ge{y) = - Gh{e* +y)- ch{e* - y) + Gh{2e*), 
where y = G'g{x). This follows from equation IjB.Sp of appendix IbI The domain of 
Ge{y) is-e* <y<e*. 

The 'Hamiltonian' form of the functional is 

{p e Alt(^,„t) I \ps\ <0:}x {s e Sym(^„t) | \ss\ < x ^ R, 
'S'l^p (P, 5, p) = {PeiPfk - P/, ) - Ge, (ps) + ssipf^ + p/, ) - Ge, {sg) 

0/ 

where the first sum is taken over all non-oriented interior edges e, and e is one of 
the oriented representatives e. The second sum is taken over all faces /. Writing it 
all out, we have 

^ifyp (P, P) - E {Mph - P/. ) + CI2 + Pe) + CI2 - ps) 

+s,{pf, + PS,) + CI2 {e: + ss) + CI2 {e: - s,) (3-10) 
-2CI2 {2ei))+Y.'^iPf. 

of 

The functional S']^p(p, s, p) is convex upwards in the variables p, s. If p and s are 
related to p by 

Pe- = G's^ [pk - Pj) = 2 arctan (tan tanh f^^-y^] ] (3.11) 



and 

.s^ = (7o ('(ii- 4- n,. "! = 2 arctan I tan I ^ 



sg = G'g^ {pk +pj)=2 arctan ( tan ( ^ ) tanh ( (Hd-fl \ \ (3.12) 
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where fj and fk are the faces to the left and right of e, respectively, then 

ShypiP^S,p) ^ S{p). 

Again, the critical points of Shyp{p) correspond to the critical points of S^yp{p, s, p) 
where p, s, and p may vary independently. (However, the domain of S^^^{p, p) 
has to be respected, of course. In particular, the constraints p e Alt(i?i„() and 
s e Sym(£'i„t) apply.) The following lemma is proved in the same way as Icmma lXTl 

Lemma 3.2. // {p,s,p) is a critical point of S^y^, then p is a critical point 
of Shyp- Conversely, suppose p is a critical point of Shyp- Define p and s by 
equations (|3.11() and (|3.12|l . Then {p,s,p) is a critical point of S^y^. At a critical 
point, Shypip) = Sl^yp{p,s,p). 

Now, we introduce the variables (p £ M,^'"*, 
instead of {p, s): 

Shyp{f,P) = S^yp(.P,s,p). 

Since 
we have 

Pe = V'e-f^e, 
Se = -<Pe~ <y5-e, 

and therefore 

SSypi'PiP) = 

X! ( ~ 2 (Ve-P/, + V'-e-P/fc ) + CI2 {91 + ^fs - ^-S) + CI2 {el - ips + ^~e) 

+ CI2 {e: - PS - p-s) + CI2 {e: + PS + <^-e) - 2 CI2 (2C)) 

of 

The domain of S-^yp{p, p) is the set of all {p, p) £ R^'"' x E,^ such that p satisfies 
the inequalities 

\ps + p-s\ < 0* and \ps - p-e\ < C- 

This is a superset of the space of hyperbolic coherent angle systems. The rela- 
tions H3.11|l and (|3.12() between p, s and p are equivalent to the relation (|2.16() be- 
tween p and p. Hence, if p and p are related by equation (|2.16(l . then S^y^{p, p) = 
'S'hyp(p)- The functional S^y^{p,p) is linear in p. As in the euclidean case, collect- 
ing the coefficients of each pf, one observes that if (/s is a hyperbolic coherent angle 
system, then does not depend on p. In fact, in that case, S'^p((/9, p) = S{p>). 

The convexity claim of thcorcm l3.1l (ii) follows from the convexity of S^yp{p, p) 
with respect to p. Equations H3.3|l are the inverse relations to equations Ij'i.lGII : see 
lemma 12.51 It is left to prove that they are compatible if (ys is a critical point of 
S under variations in the space of hyperbolic coherent angle systems. The rest 
follows from lemma 13.21 Suppose the oriented edges ei and 6*2 are both in the 
oriented boundary of a face /. We need to show that the value for pf obtained 
from the equation (|3.3() involving p±s^ is the same as the value obtained from the 
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equation involving Lp±g^ . The tangent space to the space of coherent angle systems 
is spanned by vectors of the form 

_d_ d_ 

where e\ and e2 are two oriented edges in the boundary of a face. But dS/dtpg 
is twice the right hand side of equation H3.3|l . Thus, they are consistent if is a 
critical point of ^hyp under variations in the space of coherent angle systems. This 
completes the proof of part (ii) of theorem 13. II 

3.1.4. The spherical functional. This case is similar to the hyperbolic case. 
Define the 'Lagrangian' form of the spherical circle pattern functional S'sph (defined 
in equation H2.27|) ') to be 

^,^ph : Alt(K„0 X Sym(^,„i) x ^ M, 

If, for all oriented edges e&E, 

"e- = Pf, - PS, and ws = p/, + p/, , 

where fj and fk to arc the faces to the the left and to the right of e, respectively, 
then 

Ssp\,{v,w,p) ^ S'sph (p). 
The 'Hamiltonian' form is defined as 

{p e k\i{E,nt) I be I <ei]x{se Sym(^,„t) I |SeH < ^^e} X ^ R, 
S"ph{P,s,p)^ {peiPf, - Pf,) - GeAPs) ~ ss{pf^ + Pf^) + G^-eAse) 

-Apfj + Ph)) + Y.^fPf- 

of 

Writing it all out, we have 

sii^ b, s, p) = {MPh - ) + ci2 {e: + ps) + ci2 {e: - ps) 

-Se{pf, + P/J - CI2 (0e + Ss) - CI2 {0e - Sg) 

-2CI2 {2e:)-n{pf^+pfS 



0/ 

(We have used that 012(26*) = -012(26**).) The functional S^^P^^^P) is convex 
upwards in the variables p and convex downwards in the variables s. If p and s are 
related to p by 

PS = G'e^ (pfe -Pj) = 2 arctan (^tan (^^^ tanh (^^^^-^) ) (3.14) 



ss = G'g, [pk + Pj) = 2 arctan ( tan ( ^ ) tanh ( Ph±Pl j j (3.15) 



and 

.s^ ~ n'^-. ( ni. + ~ 2 arctan I tan I , , 

where fj and are the faces to the left and right of e, respectively, then 

S"piAP,s,p) = S{p). 
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Again, the critical points of iS'sph(p) correspond to the critical points of S^Y^{p, s, p) 
where p, s, and p may vary independently. The following lemma is proved in the 
same way as lemmas |3 . 1 1 and 13.21 

Lemma 3.3. // {p,s,p) is a critical point of 3^^/^, then p is a critical point 
of Ssph- Conversely, suppose p is a critical point of Ssph- Define p and s by 
equations H3.14|l and (|3.15(l . Then {p,s,p) is a critical point of Sf^y^. At a critical 
point, Ssphip) = S^p,^{p,s,p). 

Now we introduce the variables ip G R^, 

(fs = ^{ps+ ss+tt) 

instead of {p, s): 

'S'sph(<P,P) = S"pi,ip,s,p). 

We have 

Pe= 'i^e- iP-e, 
Sg^ lfie + f-S- TI", 

and hence 

( ^ '^{VSPfj + V-SPh) + CI2 + - <P-e-) + CI2 (6** - + ip^s) 
+ CI2 - - V-e) + CI2 (C + + f-s) - 2 CI2 (20:)) 

of 

(We have used that Cl2(a;) is odd and 27r-periodic.) The domain of S^i^{ip,p) is 
the set of all {(p, p) G R^™* x 'K^ such that ip satisfies the inequalities 

\(ps+ ip-s - T^] < de and \ips ~ (p^^\ < 81. 

If ip and p are related by equation (|2.25|1 . then S^Y^{(p, p) = Ss-ph{p)- The functional 
S^^{<p, p) is linear in p. Collecting the coefficients of each pf, one observes that if 
95 is a spherical coherent angle system, then S^y^ docs not depend on p. In fact, in 

that case, S^^i-^{ip,p) = S{ip). 

Equations (|3.4|) are the inverse relations to equations (|2.25|) : see lemma 1^.101 It 
is left to prove that they are compatible if 1^9 is a critical point of S under variations 
in the space of spherical coherent angle systems. This follows in exactly the same 
way as in the hyperbolic case. The rest follows from lemma rOl This completes 
the proof of part (Hi) of theorem l3. II 

3.2. Colin de Verdiere's functionals 

Colin de Verdiere |16| considers circle packings in which the circles correspond 
to the vertices of a triangulation. He considers the 1-form 

Lu = adu + f3dv + j dw 

on the space of cuclidean triangles, where 

u=:loga;, v~logy, and w — logz 

and X, y, z and a, P, 7 are as shown in figure ITTI 
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Figure 3.1. A euclidean triangle 



It turns out that duj = 0, hence one may integrate. Define the function fao,fio.' 
on M? by 



The initial point of the integration does not matter. 

Suppose we are given a triangulation and a coherent angle system for it. Here, 
a coherent angle system is a positive function on the set of angles of the triangles, 
such that the sum in each triangle is tt, and the sum around each vertex is 2-k. For 
a function p on the vertices of the triangulation, Colin de Verdiere's functional for 
euclidean circle packings is 



where the sum is taken over all triangles, u, w, and w are the values of p on the 
vertices of each triangle, and ao, Po and 70 are the corresponding angles of the 
coherent angle system. 

The critical points of this functional correspond to the logarithmic radii of a 
circle packing. 

The hyperbolic case is treated in the same way, except that now 

u = logtanh(a::/2), w = logtanh(j//2), and w = logtanh(2/2), 

and in the definition of a coherent angle system, it is required that the sum of the 
angles in a triangle is less than tt. 

To treat circle packings with our functionals, we consider circle patterns with 
orthogonally intersecting circles on the medial decomposition of the triangulation 
(see figure in^ . 

Suppose Eq is a triangulation. Let S be the medial decomposition. Set Of, = 7r/2 
for all edges e. For simplicity, we will consider only the case of closed surfaces, so 
let $/ = 2-K for all faces /. The faces of E are of two types: those that correspond 
to faces of Sq and those that correspond to vertices of Sg . Let = Fi U F2 , where 
Fi contains the faces of the first type and F2 contains the faces of the second type. 

Consider S'^uc as function on K^i x R^^ ^nd define S2 : M^" ^ M, 



Hence, S^ip^) = Seudpi, P2), where pi G is determined as follows. Suppose 
f G Fi and fa, fb, fc G ^2 are the neighboring faces of /. Construct the euclidean 
triangle whose sides are 




(ao — a)du+ {(3o — (3) dv + (70 - 7) dw. 




S2{P2) 



min S'ouc(pi,P2)- 



eP^f. +ef^fb, e''" ^i. + e''" , and e''= /<= + e''^ . 
Let be the logarithmic radius of the inscribed circle. 



3.5. RIVIN'S FUNCTIONAL 



43 





Figure 3.2. Left: Thurston type circle pattern with overlap. Middle: Thurston 
type circle pattern with hole. Right: Add a circle to fill the hole. 



For p2 G , and pi the corresponding point in 



dpi 



-iPl,P2) = 0, 



^ (Pl,P2) = ^(P2)- 
OP2 OP2 



Consider Scdv as function on E^. It is not hard to see that ^5*2 = dScdv- This 
implies that Colin de Verdiere's euclidean functional is, up to an additive constant, 
equal to 82- 

Colin de Verdiere's functional for hyperbolic circle packings can be derived from 
Shyp in the same way. 

3.3. Digression: Thurston type circle patterns with "holes" 

A similar reduction leads to functionals for Thurston type circle patterns with 
"holes" |42| . Thurston type circle patterns were mentioned in the introduction, 
section n~T1 Here, the vertices do not correspond to intersection points. All vertices 
have three edges, but the intersection angles 9 do not have to sum to 2tt. If the sum 
is less than 2tt, the circles will overlap as in figure rOl (left). If the sum is greater 
than 27r, there will be a hole as in figure (middle). If the sum is greater than 2-k 
around all vertices, one obtains a Dclaunay type circle pattern by adding circles as 
shown in figure (right). Combinatorially, this corresponds to "truncating" all 
vertices. The intersection angles for the new edges (see figure rO| are determined 
by the equations 



'12 4- 


- ^20 4 


- ^01 


= 27r 


'23-1 


- ^30 4 


- ^02 




'31 4 


- ^10 4 


^03 


= 2ti 



so that 



923 + 6'3i), etc. 



3.4. Bragger's functional 

Like Colin de Verdiere, Bragger |12| considers circle packings in which the 
circles correspond to the vertices of a triangulation. A coherent angle system in the 
sense of section is clearly equivalent to a coherent angle system in the sense of 
this article. Thus, Bragger's functional is seen to be equal to S'/2 up to an additive 
constant. 

3.5. Rivin's functional 

Rivin |34| considers euclidean circle patterns with arbitrary prescribed inter- 
section angles. The pattern of intersection is determined by a triangulation. Circles 
correspond to faces of the triangulation. However, since Rivin allows intersection 
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Figure 3.3. A "truncated vertex". The old faces are labeled 1, 2, 3; the new face 
is labeled 0. 

angles 9 = 0, adjacent triangles may in effect belong to the same circle. Thus, 
cellular decompositions with non-triangular cases may be treated by first dissecting 
all faces into triangles. Rivin treats cone-like singularities in the vertices but not 
in the centers of the circles. His functional is up to an additive constant equal to 
5/2 with S as in equation H3.5|l . 

3.6. Leibon's functional 

Like Rivin, Leibon |27| , |28| considers circle patterns with arbitrary prescribed 
intersection angles. The circles correspond to the faces of a triangulation and 
the variables are the angles of the inscribed triangles. Whereas Rivin treats the 
euclidean case, Leibon treats the hyperbolic case. Leibon's functional cannot be 
derived directly from our functionals. In section [4. 71 wc will construct yet another 
functional, from which both Leibon's functional and S{ip) can be derived. This 
section is an exposition of Leibon's variational principle. 

Let S be a triangulation of a compact surface without boundary and let in- 
tersection angles be prescribed by a function 9 : E —t (0, tt) on the non-oriented 
edges. Assume that 9 sums to 27r around each vertex. (This assumption is made 
only to simplify the exposition. It means that we consider only the case without 
cone-like singularities. If the assumption is violated, one obtains prescribed cone- 
like singularities in the vertices. Cone-like singularities in the centers of the circles 
cannot be treated by this method.) The variables of Leibon's functional are the 
angles of the triangles of the triangulation. We associate each angle of a triangle 
with the opposite directed edge in the oriented boundary of that triangle. Thus, for 
a geodesic triangulation of a hyperbolic surface, the angles of the triangles define a 
function on the oriented edges. (In the figures, however, we label the angles in the 
usual way.) 

For a function a g E,'^ on the oriented edges, Leibon's functional is 
i/(«)=^y(al,4,4), 

where the sum is taken over all triangles i; a\,a\,a\ are the three angles of a 
triangle and 

T/(ai, a2, as) = ^ ( Cl2(2ai) + Cl2(2a2) + Cl2(2a3) 

+ Cl2(7r + ai - Q2 - 03) + Cl2(7r — ai + a2 — a^) (3.16) 
-I- Cl2(7r — ai - Q2 + Q3) + Cl2(7r — ai — a2 — as) J . 
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Figure 3.4. Left: Adjacent triangles. Right: Triangle with circumscribed circle. 
(For simplicity, euclidean triangles are drawn in the figure.) 

A coherent angle system in this setting is a function a G on the oriented 
edges, which satisfies the following conditions: 
(i) For all e e E, ag > 0. 
(a) For all triangles t ^ a\+ a\ + a\ < it. 

(in) The angles of two adjacent triangles with common edge e & E, as shown 
in figure rOI (left), satisfy 



Note that condition (Hi) implies that the angles a sum to 27r around each vertex 
(since the angles 6 do). For a Delaunay triangulation with circles intersecting 
at angles given by 9, the angles of the triangles form a coherent angle system. 
Figure rOI (right) shows that condition (Hi) is satisfied. 

Lemma 3.4 (Leibon |28| 1. The hyperbolic triangles with angles prescribed by 
a G fit together to form a Delaunay triangulation of a hyperbolic surface with 
circles intersecting at angles given by if and only if a G E.^ is a critical point of 
the functional H{a) under variations in the space of coherent angle systems. 

Leibon also shows that the functional is strictly convex upwards and uses his 
variational principle to prove an existence and uniqueness theorem for hyperbolic 
circle patterns, which is a special case of theorem II .81 128| . 

Proof of lemma 13.41 It is clear that, if the triangles may be glued together, 
one obtains a Delaunay triangulation of a hyperbolic surface with circles intersecting 
at angles given by 0. The triangles fit together if sides along which two triangles 
are supposed to be glued have the same length. Hence, the lemma follows from the 
following facts: 

(a) For each oriented edge e G E, let 

d ^ d d d 
da2 daj, do!^ dolj^ ' 

where 0,2,0^,0^2^^?, s-'"*^ figure rOI (left). The tangent space to the space of 

coherent angle systems is spanned by the vectors wg, eGE. 

(b) For positive ai, a2, 0^3 with 01+02+0:3 < tt, 
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Figure 3.5. The linear functions /, g, and g*. 



where oi is the length of the side opposite ai in a hyperbolic triangle with angles 
ai, a2, as. 

Proof of (a). Consider the linear functions 



/: 


Alt{E) - 


R 


9 •■ 




]R 


9*-- 







E 



defined by figure 1X51 (Alt(£^) denotes the space of antisymmetric functions on 
the oriented edges; see equation (|^J).) Suppose the vector spaces E,^, and M.^ 
are equipped with the canonical scalar product. Then g and g* are adjoint linear 
operators. The tangent space to the space of coherent angle systems is 

(Imr7*)^ =Kerff. 

Hence we need to show that 

Im/ = Ker5. (3.18) 

In fact, we will show that 

— > A\t{E) ^r'^ — >0 

is an exact sequence. It is easy to see that 

Ker / = 

and 

Im / C Ker g 

It is also easy to see that 

Ker .9* = 

and hence 

Im.g = (Kerg*)^ = R^ . 

Since 

dimKery = - \E\ = \E\ = dinilm/, 
equation H3.18|l follows. This proves (a). 
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Proof of (b). By equation HB.4p . 

/ d d \ , , cos °^+°^+"'^ cos 
I ^ I V{ai,a2,a3) = log ^ ^ 



V (5a2 9q;3 / ' ' sin a2 sin as 

From the hyperbolic angle cosine theorem (see, for example, Beardon [5]) 

cos a2 cos as + cos ai 

coshai = ; ; 

sma2 smaa 

one derives the half-angle formula 

„ ona "1+^2+03 „„„ -0:1+02+03 

9 di cos 9 cos o 

smh — = = . 

2 sin a2 sin 

This proves (b). ■ 

3.7. The Legendre dual of Leibon's functional 

In this section, we construct the circle pattern functional H, which is Legendre 
dual to Leibon's functional. The variables are (essentially) the edge lengths of the 
triangles. Unfortunately, we cannot give an explicit formula for H in terms of the 
dilogarithm or related functions. Whether a manageable formula exists, seems to 
be an interesting question. 

Introduce new variables P G M,^ instead of a S M,^: For each triangle, the 
angles /?2, are related to ai, a2, by 

—ai + a2 + as 





2 




ai ~ 


- a2 4 


' as 




2 




ai - 


h Q!2 - 


as 


2 



1 

P2 = "^""^"""^ (3.19) 
Note that 

_9_ _ _a_ _d_ 

dPi da2 das ' 
such that equation ()3.17|l is equivalent to 

dV 

— = 21ogsinha,. 

Let V be the Legendre transform of V, considered as a function of (3i, P2, Ps- (See 
section r^.l. in That is, let 

y ( 2 log sinh a 1 , 2 log sinh 02,2 log sinh as) = 

2/?! log sinh ai + 2/32 log sinh 02 + 2/33 log sinh 03 — V{ai,a2, as), 

where ai, a2, as and ai, 02, as are the angles and side lengths of a hyperbolic 
triangle. Then 

9l 

d (2 log sin aj ) 

For a function a £ M,^ on the non-oriented edges which satisfies the triangle in- 
equalities, define the functional 

H{2 log sinh a) = V{2 log sinh a* , 2 log sinh 02, 2 log sinh a|) — 26*2 log sinh Oe. 

teF eeE 

The first sum is taken over all triangles t, and a*, a\, a\ are the sides of triangle 
t. The second sum is taken over all non-oriented edges e. The functional H is the 
Legendre dual of Leibon's functional H: 



48 



3. OTHER VARIATIONAL PRINCIPLES 



Proposition. Suppose a g R satisfies the triangle inequalities. Then the 
hyperbolic triangles with edge lengths given by a form a Delaunay triangulation of 
a hyperbolic surface with circles intersecting at angles given by 6 if and only if 
21ogsinha is a critical point of the functional H . 

Proof. Suppose e and — e are the two orientations of the unorientcd edge 
e G E. Then the partial derivative of H with respect to 2 log sinh is 

gg(21ogsinha) 
a (2 log smh Oe) 

This partial derivative vanishes if and only if the circumscribed circles of the two 
triangles incident with e intersect at the angle Oe.] see figure (right). ■ 



CHAPTER 4 



Circle patterns and the volumes of hyperbolic 

polyhedra 

In chapter 121 we constructed circle patterns by tiling (euclidean, hyperbolic or 
spherical) 2-space with kite-shaped quadrilaterals. In this chapter, we extend these 
tilings to 'tilings' of hyperbolic 3-spacc with 3-diinensional polyhedral tiles. This 
provides a geometric interpretation of the functional S of section IXTI It turns out 
to be the sum of the volumes of the polyhedral tiles. The functional attains its 
maximum (under variation constrained to coherent angle systems) if the tiles fit 
together to form a 3-dimensional polyhedron. The maximal value is therefore the 
volume of this polyhedron. 

4.1. Schlafli's differential volume formula 

In 1852, Schlafli discovered the remarkable formula (|4.1|l for the differential of 
the volume of a spherical simplex of arbitrary dimension |38| . In 1936, Kneser gave 
an elegant proof (in a nasty journal) |25) . It works also for hyperbolic simplices. 

Theorem 4.1 (Schlafli's differential volume formula). For an n-dimensional 
simplex S in spherical or hyperbolic space, let V{S) be its n-dimensional volume. 
Denote the {n — l)-faces by Si, < i < n. Let Ujk be the {n — 2) -dimensional 
volume of the {n — 2)- face Sj H Sk, and let ajk be the the dihedral angle between Sj 
and Sk- Then the differential of the volume function V is 

dV ^ —^-^ ^ Ojkdajk, (4.1) 

0<j<k<n 

where e = I in the spherical and e = — 1 in the hyperbolic case. 

For 3-dimensional hyperbolic simplices (tetrahedra), formula (|4.1|) may be writ- 
ten 

dV^~l J2 (4-2) 

edges j 

where the sum is taken over the edges j, and aj, aj are the length and dihedral 
angle of edge j. In fact, this formula holds for arbitrary polyhedral shapes: 

Corollary. Formula H4.2|l holds for arbitrary compact hyperbolic 3-manifolds 
with polyhedral boundary (in particular, for 3-dimensional hyperbolic polyhedra), 
when they are deformed without changing the combinatorial type of the polyhedral 
boundary. 

This follows by triangulating the 3-manifold with polyhedral boundary. At 
all edges of the triangulation, except those which are contained in edges of the 
polyhedral boundary, the sum of dihedral angles is constant. 

Now consider hyperbolic polyhedra (or manifolds with polyhedral boundary) 
with some vertices in the infinite boundary of H^. Even though they are not 
compact, their volume is finite. However, formula H4.2|l does not apply directly. 
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because some edges have infinite length. The following lemma is ascribed to Milnor 
m P- 576]. 1 

Lemma 4.1 (Milnor). Suppose P is a S- dimensional hyperbolic polyhedron with 
some vertices on the infinite boundary. For each infinite vertex of P , choose an 
arbitrary horosphere centered at this vertex. For a finite edge j of P , let Oj be its 
length. For an edge j with one ( or two ) vertices at infinity, let Oj be its length up 
to the cut-off point (s) at the corresponding horosphere (s). (This may be negative if 
the horospheres are not small enough.) With this agreement, the following holds: 
If P is deformed in such a way that its combinatorial type does not change and the 
infinite vertices stay infinite, then the volume differential is given by formula H4.2|l . 

Proof. Because the sum of dihedral angles at an infinite vertex with n edges is 
constantly {n — 2)tt during the deformation, the right hand side of equation (|4.2() is 
independent of the choice of horospheres. We may thus assume that the horospheres 
are small enough, so that each intersects only the planes incident with the vertex 
at which it is centered and no other planes or horospheres. By triangulating the 
polyhedron appropriately, one finds that it suffices to consider a triply orthogonal 
tetrahedron with one vertex at infinity, see figure [HTI (left). The volume of such a 
tetrahedron can be found by a straightforward integration without using Schlafli's 
formula, as shown in appendix [0 Differentiate equation l)C.l|l to obtain, after a 
straightforward calculation. 



dV 1, / cos2/3 

— = --logWl 

oa 2 V cos^ a 

Truncate the tetrahedron by the horosphere centered at the infinite vertex which 

touches the opposite face. Then the truncated sides with dihedral angles a and 

J — a have length a = and a = — log y^l — f^fs-^, respectively. Hence, 
dV\ ^ — — (ada + adiiT ~ a)). 

I l3=const. 2 ^ ' 

This proves the theorem. ■ 
4.2. A prototypical variational principle and its Legendre dual 

In this section, we take a little detour to illustrate the fundamental ideas which 
are applied in the next sections. It may be skipped, or read after the following 
sections. We derive a pair of Legendre dual variational principles connected with 
triangulations of hyperbolic 3-manifolds with polyhedral boundary. The arguments 
below may be adapted to deal with the case where all vertices of the triangulation 
are at infinity |39| . In that case, one does obtain convex functionals. 

Let T be a finite triangulation of a compact topological 3-manifold M with 
boundary, and let To, . . . , Ta be the sets of vertices, edges, triangles, and tetrahcdra. 
Let S be the manifold of shapes of hyperbolic tetrahcdra. This is diffeomorphic to 
a connected open subset of R^. Global coordinates are the six dihedral angles (or, 
alternatively, the six edge lengths), which satisfy certain inequalities. 

For a function tj : — > 5, which assigns a shape to each combinatorial tetra- 
hedron of T independently, let 

'S'proto(f) ~ 2^y(a(i)), (4.3) 

tela 

where V is the volume function on S. This defines a differentiable functional 
"Sproto ■ S ^ ^ M. 



"'^Howevcr, the proof given by Rivin does not seem to be valid, because, in the "slicing a 
carrot" argument, Schlafli's formula is applied to shapes bounded by planes and horospheres. 



4.3. THE EUCLIDEAN FUNCTIONAL 
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Let $ : Ti n dM — > ]R+ be a function that assigns a positive angle to each 
boundary edge. A shape assignment ct : 73 — > 5 is coherent, if the dihedral angles 
sum to 27r around each interior edge and to <&(e) around each boundary edge e. 

Proposition. A coherent shape assignment u G S'^^ is a critical point of the 
functional Sproto under variations in the space of coherent shape assignments if and 
only if the hyperbolic tetrahedra 17(73) fi-t together to form a hyperbolic manifold with 
polyhedral boundary, where the boundary angles are prescribed by $. 

Proof. Let a- be a coherent shape assignment. Suppose ii,t2 G share an 
edge e € Ti. Let ai, 02 be the dihedral angles of (T(ii), cr(i2) at e, and let ai, 02 be 
the side lengths of e in a(ti),a{t2). There is a variation v in the space of coherent 
shape assignments with 1 = dai{v) = —da2{v) which keeps all other dihedral angles 
constant. By Schlafli's formula (|4.2(l . 

dSp,oto{v) = -ai + 02- 

Hence, if cr is a critical point of S'proto under variations in the space of coherent 
shape assignments, then cr assigns a unique length to all edges e G 7i. Therefore, 
the tetrahedra fit together. 

The converse statement follows, because variations like v span the tangent space 
to the space of coherent shape assignments. I 

To define the Legendre dual functional, consider a different subset of shape 
assignments: the space TZ C S^-^ of those a £ S^^ which assign a unique length to 
each edge e € Ti. This means, if ti, i2 G Tij are both incident with e, then the length 
of e in a(ti) equals the length of e in a{t2)- For each such shape assignment a Cz TZ 
and each edge e e 7i, let ae(cr) be the length assigned to e by a. The functions 
Ce : 7?. — > R-i- (which satisfy certain inequalities) form a system of coordinates of TZ. 
Hence, TZ is diffeomorphic to an open subset of . For a tetrahedron t £7^ with 
edge e S t n 71, let at.e{<^) be the dihedral angle at e in (j{t). Extend the function 
which was defined on the boundary edges, to a function on all edges by setting 
$(e) = 27r if e is an interior edge. 

Define the functional S'proto '■ TZ —^ 

Voto(CT) ^-Y.{ ae{cT)at,eia) + 2F(a(0)) + *(e)ae(^T). 

teTs eetnTi eSTi 

Proposition. A shape assignment G TZ is a critical point of the functional 
Sproto (under variations in TZ) if and only if the tetrahedra (7(73) fit together to 
form a hyperbolic manifold with polyhedral boundary, where the boundary angles 
are prescribed by $. 

Proof. By Schlafli's differential volume formula H4.2|l . 

dSproto = ^ (*(e) - ^ at,e) doe. 

t3e 

Thus, at a critical point, the sum of all dihedral angles at an edge e is $(e). ■ 

4.3. The euclidean functional 

Suppose that < 9 < tt and ipi > 0, (p2 > 0. If 

(fil + IP2 = - 0, 

we define a hyperbolic polyhedron Pg{ipi, ip2) as follows. 

If 1^1 and ip2 are not greater than ^, let Pe(vii ^2) be the hyperbolic polyhedron 
shown in figure 14.11 in the Poincare half-space model of hyperbolic 3-space: = 
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Figure 4.1. The hyperbolic polyhedron Pb[lpi,lp2) for ipi + Lp2 = 0* (where 6* = 
TV — 9) shown in the Poincare half-space model. The dihedral angle at the four 
unmarked edges is ^. 




Figure 4.2. Left: The polyhedron of figure Wa\ projected orthogonally into the 
z-plane. Right: The corresponding picture for the case ip2> 



{{x,y,z) e H"^ I z > 0} with metric \ds\ = \\/ x'^ + y'^ + z^. It has three infinite 
vertices: two lie in the plane z = and third is the infinite point compactifying the 
closed half-space. The dihedral angles are as indicated. The unmarked edges have 
dihedral angle \. 

If either ip\ or ipi are greater than ^, we define Pe{}p\, ip2) as an algebraic sum. 
Figure IT!^ (left) shows the orthogonal projection of figure EtI into the z-plane. The 
corresponding figure for the case (/?2 > § is shown on the right. In this case, let 
Pi be the polyhedron bounded by the planes corresponding to the circle around A 
and the line segments AB, BD, and DA; and let P2 be the polyhedron bounded 
by the planes corresponding to the circle around C and the line segments CB, BD, 
and DC. Define Pg{ifii, (^2) as the algebraic sum Pe{ipi, (^2) = Pi — P2- Extend the 
volume function linearly to such algebraic sums: V{Pi — P2) = ^(^'i) — V{P2)- 
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Proposition. Suppose G E,^'"* is a euclidean coherent angle system (see 
section \2. Then the functional S{(p) of theorem \S.l\ satisfies 

where V is the volume function, the sum is taken over all non- oriented interior 
edges e, and e, — e are the two corresponding oriented edges. 

Proof. Because of equation H3.5(l for S{if), where is a euclidean coherent 
angle system, the proposition follows directly from the following lemma. ■ 

Lemma 4.2. The volume of Pg{Lpi,ip2) is 

V{Pe{ipu<P2)) - \ Ch{2ipi) + i Cl2(2</J2). (4.4) 

Proof. The polyhedron Pg{ipi,ip2) may be decomposed into four triply or- 
thogonal tetrahedra (see appendix [n]) , each having two vertices at infinity. The 
characteristic angle is tpi for two of the tetrahedra and ip2 for the other two. Equa- 
tion (|4.4() follows from equation (jC.2p . (This argument works also in the case where 
(ySi or iy92 arc greater than j.) M 

Corollary. If ip G M.^'"* is a critical point of S under variations in the space 
of euclidean coherent angle systems, then the polyhedra Pg^ {ipg, (f-g) fit together 
to form a hyperbolic orbifold M with polyhedral boundary. M is a manifold with 
polyhedral boundary if $y = 27r for all interior faces f £ F. The volume of M is 

V{M) = ^(^S{v)+ ^ Cl2(20:)) =i(5e.e(p)+ ^ Ch{20:)), 
where p G is the corresponding critical point of Seuc- 

4.4. The spherical functional 

In the previous section, we defined Pg(tpi, ip2) for pi + ip2 ^ tt — 0. If 

n — 6<'p>i+(p2<TT + 9 and | (/'i — <y92 1 < tt — 6, 

let Pg{ipi, 1P2) be the hyperbolic polyhedron shown in figure^31in the Poincarc ball 
model of hyperbolic space: 

{ix,y,z) G I + _^ < 1} 

with metric 

2 

1 — a;-^ — y-^ 

(If ipi or ip2 are greater than ^, then Pg{ipi, P2) is defined by an algebraic sum, as 
in the previous section. We will not dwell on this point.) 



\ds\ = \/ dx"^ + dy"^ + dz"^ 



Proposition. Suppose p) G R^'"' is a spherical coherent angle system (see 
section \2. ?| ). Then the functional S{p)) of theorem VS . 1\ satisfies 

S{^) = ^V{PgS^,,^^,)) 

where V is the volume function, the sum is taken over all non-oriented interior 
edges e, and e, —e are the two corresponding oriented edges. 
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Figure 4.3. The hyperbolic polyhedron Pg{Lpi,(p2) for n — 9 < Lpi + ip2 < n + 
and |<pi — </'2| < TT — 9, shown in the Poincare ball model. The dihedral angles at 
the edges marked ai and a2 are 2ipi and 2ip2. The other two edges of the vertex 
marked 'o' (which is at {x,y,z) = (0,0,0)) have angle 9; and the faces opposite 
vertex 'o' intersect at the angle 9* — it — 9. The four remaining edges have dihedral 
angle ^ . The polyhedron has two vertices in the infinite boundary. 

Corollary. If (p E M.^'"' is a critical point of S under variations in the space 
of spherical coherent angle systems, then the polyhedra Pg^ [ipg, (p-e) fi-t together 
to form a hyperbolic orbifold M with polyhedral boundary. M is a manifold with 
polyhedral boundary if = 27r for all interior faces f £ F. The volume of M is 

V{M) = i5(^) = ^Ssphip), 

where p £ is the corresponding critical point of Sgph- 

The proposition follows directly from the following lemma. 

Lemma 4.3. The volume of Pg{pi, P2) is 

V{Peipi,P2)) = I ( Cl2(r + pi- P2) + Cl2{9* -pi+ P2) 

2V (4.5) 

+ Cl2(r +PI+ P2) + C\2{9* -pi-p2)-2 Cl2(20*) j . 

Proof of Lemma ^31 We will show that the derivatives of both sides of 
equation H4.5() with respect to pi and p2 are equal. Since both sides tend to 
as (^1 ^ and p2 9*, equation H4.5|l follows. 

Consider deformations of Pg{pi, P2) during which all dihedral angles except 
2pi and 2p2 remain constant. By Milnor's generalization of Schlafli's differential 
volume formula flcmma l4.1() . 

dV{Pg{pi,p2)) = -aidpi -a2dp2, 

where ai, a2 are the lengths of the sides with dihedral angles 2pi and 2(p2- (And 
this is also true in the case where Pe{(pi,p2) is defined as an algebraic sum.) By 
equation ljB.4|l . the differential of the right hand side of equation ()4.5|l is 

pi dpi + p2 d(p2, 

where pi and p2 are given by equation (|3.4(l . Let ri and r2 be the spherical radii 
of the circles in the infinite boundary shown in figure 14.31 Then ri and r2 are 
related to pi and p2 by equation H2.23|l . This implies pi = — ai and p2 = —02; see 
figure ^31 This completes the proof. ■ 
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Figure 4.4. In a hyperbolic right-angled triangle with one vertex at infinity, the 
length a and the angle r satisfy the equation a = — log tan ^. 




Figure 4.5. The hyperbolic polyhedron Pg{Lpi, (^2) for (fii + (p2 < tt — 9, shown in 
the Poincare ball model. The dihedral angles are as indicated, where 6* ~ n — 0. 
The unmarked edges have dihedral angle ^. The polyhedron has four vertices in 
the infinite boundary. It is symmetric with respect to reflection on the equatorial 
plane. 

4.5. The hyperbolic functional 

Now we define the hyperbolic polyhedron Pg(ipi, ip2) for 

Lpi+(p2<TT — 9. 

In this case, let Pg{ipi, ip2) be the hyperbolic polyhedron shown in figure 1^751 in the 
Poincare ball model of hyperbolic space. (If tpi or ip2 are greater than ^, then 
Pe{ipi,ip2) is defined by an algebraic sum, as described in section ing We will not 
dwell on this point.) 

Proposition. Suppose ip £ R^™' is a hyperbolic coherent angle system (see 
section \2. 7| ). Then the functional S(if) of theorem \S.l\ satisftes 

where V is the volume function, the sum is taken over all non- oriented interior 
edges e, and e, —e are the two corresponding oriented edges. 

Corollary. If (p E M,^'"' is a critical point of S under variations in the space 
of hyperbolic coherent angle systems, then the polyhedra Pg^ {ipg, ips) fit together 
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to form a hyperbolic orbifold M with polyhedral boundary. M is a manifold with 
polyhedral boundary if = 27r for all interior faces f G F. The volume of M is 

where p G is the corresponding critical point of Shyp ■ 

The proposition follows directly from the following lemma. 
Lemma 4.4. The volume of P0{ipi, tp2) is 

v{Peiipuip2)) =ci2(r + </52) + ci2(r-v^i + (^2) 

+ ci2(r + ipi + ip2) + ci2(r - - ^2) - 2 ci2(2r). ^ ■ ' 

Lemma 14.41 is proved in the same way as lemma ^31 

Proof of Lemma ^3 We will show that the derivatives of both sides of 
equation H4.6|) with respect to ipi and are equal. Since both sides tend to 
as (/?! — > and Lp2 —^6*, equation H4.5|l follows. 

Under deformations of Pg{ipi, ip2) during which all dihedral angles except 2tpi 
and 2ip2 remain constant, the volume differential is 

dV{Pei(pi,(p2)) = -aid(pi ~ a2dip2, 

where ai, 02 are the lengths of the sides with dihedral angles 2(pi and 2(^2- (And 
this is also true in the case where Pe{(pi,ip2) is defined as an algebraic sum.) By 
equation ljB.4p . the differential of the right hand side of equation (|4.t)|) is 

2pi difi + 2p2 dip2, 

where pi and p2 are given by equation The symmetry plane intersects 

Pe{ipi,ip2) in a hyperbolic kite-shaped quadrilateral with angles 9, (pi, tp2- Let 
the hyperbolic side lengths of this quadrilateral be ri and r2 . Then ri and r2 are 
related to pi and p2 by equation H2.14|l (see lemma 1^5]) . Furthermore, Oj ~ —2pj 
(see figure ^3J- This completes the proof. ■ 

4.6. Leibon's functional 

Schlafli's differential volume formula is the basis for Leibon's functional (see 
section r3.6|l as well. Lemma |4.5I and Schlafli's differential volume formula provide 
another proof of Leibon's variational principle, lemma 

An ideal triangular prism is a convex hyperbolic polyhedron with six ideal 
vertices, which is symmetric with respect to a hyperbolic plane (see figure ETjl . An 



4.6. LEIBON'S FUNCTIONAL 



57 




Figure 4.7. An ideal triangular prism, shown in the Poincare ball model. The 
equatorial plane is the plane of symmetry. 





72 log sinh | V 











Figure 4.8. Truncated quadrilateral face of the ideal prism. 



ideal triangular prism intersects the plane of symmetry in a hyperbolic triangle. 
The angles ai, a2, of this triangle determine the shape of the prism. They 
are also the dihedral angles of the three edges intersecting the symmetry plane. 
The other dihedral angles are ^ — (3i, ^ — (32, and ^ — (S^, with the Pj defined by 
equations (|3.19|) . 

Lemma 4.5 (Leibon |27j)28| ). The hyperbolic volume of the ideal prism shown 
in figure \4-. ?| is V{ai,a2,ci'i) as defined by equation H3.16|l . 

Truncate the ideal prism by horospheres which are centered at the ideal vertices 
and touch the symmetry plane. The truncated edge lengths are for the edges which 
intersect the symmetry plane, and 2 log sinh ^ for the other edges, where aj is the 
length of the corresponding side of the hyperbolic triangle with angles a\, 012, 0^3 ■ 
(See figure \J~^ ) 
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Figure 4.9. Left: Schematic picture of one of the three pieces into which the ideal 
prism of figure UTI is decomposed. It has four ideal vertices, which are marked 'o'. 
The dihedral angles at four of the edges are ^ . The other dihedral angles are as in- 
dicated. Right: The triangular intersection with the symmetry plane is decomposed 
into three isosceles triangles. 



4.7. A common ancestor of Leibon's and our functionals 

The geometric interpretation of both Leibon's functional and the functional 
S(ip) in terms of hyperbolic volume clarifies their mutual relationship. In this sec- 
tion, we construct another functional, T{ip,P), from which both Leibon's functional 
and S{ip) can be derived. Let S be a cell decomposition of a compact surface with- 
out boundary, let intersection angles be prescribed by a function 9 : E ^ {0,n) on 
the non-oriented edges, and let cone angles in the centers of the circles be prescribed 
by a function <I> : F ^ (0, 27r) on the faces. 

The ideal prism shown in figure ETI can be decomposed into three polyhedral 
pieces like the one shown in figure The hyperbolic polyhedron Pg{(pi, (^2) (see 
figure can be decomposed into two pieces of this kind. Let W{(f,l3) be the 
hyperbolic volume of such a polyhedron. 

A coherent angle system in this setting is a pair (p E R^, /3 G M.^ of functions 
on the oriented edges satisfying 

(i) Foralleei?, 

<fs > 0, I3g > 0, and ipg + (3g<^. (4.7) 

(ii) For all faces f € F, 

Y.2^s^^f, (4.8) 

where the sum is taken over all oriented edges in the oriented boundary of /. 
(Hi) For all e£E, 

l3s+(3-e^9,. (4.9) 
For (f G R^, P G which satisfy (i), define the functional 

(Condition (i) has to be fulfilled for the corresponding polyhedron of figure 
(left) to exist.) 

Proposition. A coherent angle system ip, P is a critical point ofT(p, /3) under 
variations in the space of coherent angle systems, if and only if the isosceles triangles 
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with angles 2ipg and Pg as in figure (right) fit together to form a circle pattern 
with intersection angles 9 and cone angles $ in the centers of circles. 

Proof. The tangent space to the space of of coherent angle systems is spanned 
by vectors of the form 

d d 



dPe 9/9-, 

and 

d d 



where ei and 6*2 are two consecutive oriented edges in the oriented boundary of the 
same face. If and only if 

then the polyhedral pieces corresponding to the edges e and — e fit together to form 
a polyhedron Pg^{ifg^,ifg^) (see figure 
If and only if 

^ ^ V(^,/3) = 



for all pairs of consecutive edges Cj, Ck in the oriented boundary of a face, then the 
corresponding polyhedral pieces fit together to form a (not necessarily triangular) 
ideal prism as shown in figure ^21 (with a cone- like singular line if 7^ 27r). ■ 

The following two propositions explain how the functional S{ip) and Leibon's 
fimctional can be obtained as two different reductions of T{ip,P). 

Proposition. Ifip, (3 is a coherent angle system in the sense of this section, 
then ip is a hyperbolic coherent angle system in the sense of section \2. 7| 

Suppose If G is a hyperbolic coherent angle system in the sense of sec- 
tion \2. 7j Then there exists exactly one (3{(p) £ M.^ which satisfies (|4.7|l and H4.9|l . 
such that Lp, f3{ip) is a critical point of T{ip^ f3) under variations with fixed ip, and 
with j3 constrained by equation H4.9|l . The hyperbolic circle pattern functional is 

S{p)=T{p>,(i{p,)). 

(Sketch of a) proof. The first claim about coherent angle systems follows 
directly from the definitions. 

Now, suppose ip £ M.^ is a hyperbolic coherent angle system in the sense of 
section For each non-oriented edge e = {e, — e}, construct the hyperbolic kite 
with angles 2ipg, 2ip-g, and 9^, as in figure (left). The angle 9^ is divided by one 
of the diagonals into two angles (3g and (3-g, such that (p, f3 form a coherent angle 
system in the sense of this section. This {(p, /?) is the unique critical point of T{(p, (3) 
under variations with fixed ip, and with /3 constrained by equations H4.9|l . Indeed, 
it follows from Schlafli's formula, that for every such critical point, the polyhedra of 
figure IT^ (left) fit together in pairs to form polyhedra as in figure IT^ (This means 
that the isosceles triangles of figure f right ) fit together to form kites.) Equation 
S{ip) = T{tp, P{ip)) follows form the additivity of the volume function. ■ 

Proposition. Suppose that S is a triangulation, 9 sums to 27r around each 
vertex, and <I>j = 27r for all faces f . 

If p, (3 is a coherent angle system in the sense of this section, then a G M.^ , 
defined by equations H3.19(l . is a coherent angle system in the sense of section XS.fA 

Suppose a G E,^ is a coherent angle system in the sense of section \8.f)\ and 
(3 G is defined by equations (|3.19() . Then there exists exactly one (p{(3) G 
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which satisfies 1)4. 7() and (|4.8|l . such thatip, is a critical point ofT((p,l3) under 

variations with fixed (3, and with Lp constrained by equation H4.8|l . Leibon's circle 
pattern functional is 

H{f3)=Tiipip),f3). 
The proof is similar to the proof of the last proposition. 
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A computer implementation 

The variational method of constructing circle patterns described in chapter |21 
is well suited for computer implementation. Here, we present an implementation in 
Java. The classes in the cellularSurf ace and circlePattern package hierarchies 
were written by the author. The mfc, numericalMethods, and moebiusViewer 
package hierarchies are part of the jtem project j.241 . Of these classes, the author 
has written ChebyshevApproximation, Clausen, and HermitianCircle. 

The render package contains Ken Perlin's 3D renderer |32| . Originally, the 
author used Oorange, an experimental programming tool for Java and a part of 
the jtem project, to experiment with the circle pattern classes. Here, we use the 
BeanShell [4], a Java scripting tool, to illustrate how the classes work. 

5.1. Getting started 

Unix. Download the zip-file circlepatternsoftware.zip from the URL 
: / / www ■ math . tu-berlin . de/ ~ springb/ software! ^^'^ unzip it. Change into 
the directory circlepatternsof tware and execute the shell script bin/bsh with 
one of the example BeanShell scripts in the examples directory as parameter: 

> bin/bsh exsmiples/cube . bsh 

You should see some text output like this 

Surface has 6 faces, 12 edges, and 8 vertices. 
Minimizing ... Minimum found after 237ms. 
Laying out circles ... Done. 

and one or more windows displaying circle patterns should open (see figure IFTTl . 

Windows. Download the file circlepatternsoftware.zip and unzip it as 
above. Open the circlepatternsof tware folder and, in it, the folders bin and 
excimples. Drag and drop an example BeanShell script (such as cube.bsh.) from 
the examples folder onto the bsh.bat batch file in the bin folder. 

Alternatively, you may start a command shell and invoke the batch file bsh . bat 
with one of the example BeanShell scripts as argument. 

What the scripts do. The scripts bsh and bsh.bat only call the BeanShell 
interpreter with the right class path, passing on the argument. Should they not 
work properly — for example, because you have a Mac — here is what to do. Make 
sure that the directory els and the jar-file(s) in the directory lib are in the class 
path. (Of course, no other versions of the BeanShell, Ken Perlin's renderer, or the 
jtem classes should be in the class path.) Invoke a Java virtual machine and run 
the class bsh. Interpreter with a BeanShell script as argument. 

Running the BeanShell in interactive mode. If you execute the scripts 
bsh or bsh . bat without arguments, the BeanShell will run in interactive mode. You 
can run BeanShell scripts using the BeanShell command source () . Listing^shows 
an example session with the BeanShell in interactive mode. After the BeanShell is 
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invoked, it prompts the user to enter commands. In the example session, the script 
cube.bsh is run and then the gradient of the functional is printed out. 

5.2. The example scripts 

Figure shows circle patterns produced by BeanShcU scripts in the exsunples 
directory. The scripts cube.bsh, dodecahedron. bsh and icosahedron.bsh pro- 
duce the obvious trivial circle patterns in the sphere. A slightly more interesting 
pattern is produced by truncated_cube .bsh. All circles intersect at right angles. 
Half the circles form a packing with the combinatorics of a truncated cube. The 
script branched. bsh produces two circle patterns. The first is the orthogonal pat- 
tern consisting of two packings with icosahedral and dodecahedral combinatorics. 
The second is a branched pattern with the same combinatorics and intersection 
angles. There are two opposite branched faces with cone angle An. 

We will illustrate how the scripts work by taking a closer look at cube . bsh, see 
listing[5| First, the relevant classes are imported. 

1 import cellularSurf ace . *; 

2 import cellularSurf ace . examples . *; 

3 import circlePattern. variational.*; 

4 import circlePattern .viewer .perlin. *; 

5 import circlePattern. viewer .moebius.*; 

The package cellularSurf ace contains only one class, CellularSurf ace, which 
implements a combinatorial model for cell decompositions of surfaces. The pack- 
age cellularSurf ace . examples contains classes like Cube, which inherit from 
CellularSurf ace but have different constructors. Thus, 

7 CellularSurf ace surface = new Cube(); 

assigns to the variable surface a CellularSurf ace representing the combinatorial 
type of a cube. 

The package circlePattern . variational contains the classes SphericalData 
and SphericalLayout. The class SphericalData holds the intersection angles, 
cone angles, and variables p for a spherical circle pattern. It can compute the 
correct p by minimizing the spherical circle pattern functional S'sph(p). Given an 
instance of SphericalData, the class SphericalLayout calculates the positions 
of the centers and intersection points of a spherical circle pattern. (To construct 
euclidean or hyperbolic circle patterns, use the analogous classes EuclideanData, 
EuclideanLayout, HyperbolicData, and HyperbolicLayout.) 

8 GenericData data ~ new SphericalData(); 

9 GenericLayout layout — new SphericalLayout(); 

The next lines print out some information about the surface. 

11 System. out. println(" Surface has " -I- surf ace. getNumFaces() 

12 + " faces, " -I- surface. getNumEdges() 

13 -I- " edges, and " + surf ace. getNum Vertices () 

14 -I- " vertices. " ); 

Next, we set up the data. Since three edges meet at each vertex, the exterior 
intersection angles 9 are set to |7r. All $ are set to 2tt because there are no 
boundary faces and the pattern is unbranched. 

16 data.setSurface ( swr/ace ); 

17 data.assignAllCapitalPhi(2.0*Math.PI); 

18 data.assignAllTheta(2.0/3.0*Math.PI); 
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> bin/bsh 

BeanShell 1.2.7 - by Pat Niemeyer (pat9pat.net) 

bsh 7o source("examples/cube.bsh"); 

Surface has 6 faces, 12 edges, and 8 vertices. 

Minimizing ... Minimum found after 232ms. 

Laying out circles ... Done. 

bsh 7. for (int i = 0; i < 6; i++) { 

System. out. println(data.getGradient(i)); 

} 

4 . 5469956688748425E-8 
1 . 9296370279420216E-8 
1 . 9296370279420216E-8 
1 . 9296370279420216E-8 
1 . 9296370279420216E-8 
4 . 5469954468302376E-8 
bsh 7o exit(); 
> 

Listing 1. Example session: Running the BeanShell interactively. 



1 import cellularSurf ace . *; 

2 import cellularSurf ace . examples . *; 

3 import circlePattern. variational.*; 

4 import circlePattern. viewer .perlin. *; 

5 import circlePattern. viewer .moebius.*; 

6 

7 CellularSurf ace surface = new Cube(); 

8 GenericData data = new SphericalData(); 

9 GenericLayout layout = new SphericalLayout(); 

10 

11 System. out. println(" Surface has " + surf ace. getNumFaces() 

12 + " faces, " + surface. getNumEdges() 

13 + " edges, and " + surf ace. gctNum Vertices () 

14 + " vertices. " ); 

15 

16 data.setSurface ( surface ); 

17 data.assignAllCapitalPhi(2.0*Math.PI); 

18 data.assignAllTheta(2.0/3.0*Math.PI); 

19 

20 //data. setRho (0, 1.0); 

21 

22 System. out. print(" Minimizing ... " ); 

23 data.adjustRho(); 

24 

25 System. out. print("Laying out circles ... " ); 

26 layout. layout ( data); 

27 System.out.println("Done."); 

28 

29 PerlinView.show(/a?/OMt); 

30 //Moebius View. show(layout); 

Listing 2. cube. bsh 
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cube.bsh dodecahedron . bsh icosahedron.bsh 




truncated cube.bsh branched. bsh 




hex_grid . bsh 




quadmesh.bsh hyperbolic_mesh.bsh 



Figure 5.1. Circle patterns produced by the example BcanShell scripts. 
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Then, data calculates the correct p, and from this, layout calculates the coordmates 
of the centers and intersection points. 

22 System. out. print("Minimizing ... " ); 

23 data.adjustRho(); 

24 

25 System. out. print("Laying out circles ... " ); 

26 layout. layout ( data) ; 

27 System.out.println("Done."); 

Finally, wc display the resulting pattern in a 3-dimcnsional spherical view. 

29 PcrlinVicw.show(Za?/oMt); 

Dragging the mouse rotates the sphere. 
If you uncomment the last line, 

30 //MocbiusVicw.show(layout); 

another window will pop up with a view of the same circle pattern projected stereo- 
graphically to the plane. Dragging the mouse will translate the pattern. Shift click 
inside the window to select between translate, rotate, scale, info-coord, and three- 
point-transform mode. In three-point-transform mode, you can Mobius transform 
the image by dragging three points. Press the i-kcy to perform an inversion on the 
unit circle. 

The resulting circle pattern is fairly symmetric — all radii arc about equal — 
even though the solution is only unique up to Mobius transformations. This is so 
because initially, before the spherical functional is minimized, all p are equal to 0. 
If you uncomment line 20, 

20 //data. sctRho (0, 1.0); 

you will get a more unsymmetric solution. 

5.3. Class overview 

Below, we list the packages and most important classes with some comments. 
For more detailed documentation, sec sections 15 .41 and 15. 51 and the code documen- 
tation in the javadocs directory. 

classes pertaining to cellular surfaces 

package cellularSurf ace : contains only one class. 

class CellularSurf ace : implements a winged edge model for cell decomposi- 
tions of surfaces. 

package cellularSurf ace . examples : contains classes which inherit from the 
class CellularSurf ace. Their constructors create specific cell decomposi- 
tions. 

class Cube extends CellularSurf ace : the combinatorial type of the cube, 
class Dodecahedron extends CellularSurf ace : the combinatorial type of 
the dodecahedron. 

class ProjectivizedCube extends CellularSurf ace : the combinatorial type 
of a cube with diametrically opposite points identified. A non-orientable cell 
decomposition with 3 faces, 6 edges, and 4 vertices. 

class ProjectivizedDodecahedron extends CellularSurf ace : the combi- 
natorial type of a dodecahedron with diametrically opposite points identi- 
fied. A non-orientable cell decomposition with 6 faces, 15 edges, and 10 
vertices. 
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class QuadMesh extends CellularSurf ace : a cell decomposition of the disc 
consisting of a rectangular grid of quadrilateral faces. (No boundary edges.) 

class HexGrid extends CellularSurf ace : a cell decomposition of the disc 
consisting of a triangular grid of hexagonal faces. (No boundary edges.) 

classes pertaining to circle patterns 

package circlePattern. variational : 

class CPMath : provides static methods for mathematical functions that are 
used in different classes. 

class EuclideanFunctional : provides a static method to compute the value of 
the euclidean functional S'ouc(p) and its gradient, given a CellularSurf ace, 
intersection angles, cone angles, and p. 

class HyperbolicFunctional : same for the hyperbolic functional S'hyp(p). 

class SphericalFunctional : same for the spherical functional S'sph(p). 

abstract class GenericData: abstract superclass of the following three classes. 
Properties: a CellularSurf ace, intersection angles, and cone angles. 

class EuclideanData extends GenericData : minimizes the euclidean func- 
tional. 

class HyperbolicData extends GenericData: minimizes the hyperbolic func- 
tional. 

class SphericalData extends GenericData : minimizes the spherical func- 
tional. 

abstract class GenericLayout : abstract superclass of the following three 
classes. 

class EuclideanLayout extends GenericLayout : responsible for laying out 

the circles, given an instance of EuclideanData with correct p. 
class HyperbolicLayout extends GenericLayout: same for HyperbolicData. 
class SphericalLayout extends GenericLayout : same for SphericalData. 
package circlePattern. viewer .moebius : helper classes to view circle patterns 

in the MoebiusViewer. 
class CirclePatternShape : This is an adapter class to view circle patterns in 

the MoebiusViewer. Holds a reference to a GenericLayout and implements 

the interface MoebiusShape. 
class MoebiusView : the static method sliow(GenericLayout layout) opens a 

frame with a MoebiusViewer displaying the circle pattern layout. 

package circlePattern. viewer .perlin : contains only one class. 

class PerlinView : The static method show(GenericLayout layout) opens a 
frame showing a 3-dimcnsional spherical view of the circle pattern layout in 
Ken Pcrlin's renderer. 

mathematical foundation classes 

package mfc. number : provides the class Complex for complex numbers, 
package mfc. vector : provides the classes Real3 and Complex2 for real 3-vectors 
and complex 2- vectors. 

package mfc. matrix : provides classes for different types of complex matrices, 
abstract class AbstractComplex2By2 : abstract superclass of the other matrix 

classes. Provides functionality, but access is restricted, 
class Complex2By2 extends AbstractComplex2By2 : class for general complex 

2x2 matrices. 

class HermitianComplex2By2 extends AbstractComplex2By2 : class for Her- 
mitian complex 2x2 matrices. 

package mfc. group : provides classes for matrix groups. 
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class Moebius extends AbstractComplex2By2 : class for elements of the group 
of orientation preserving eonformal self-maps of the Riemann sphere. 

package mfc . geometry : provides classes for some geometric objects. 

class ComplexProjectivel extends Complex2 : homogeneous coordinates of 

1-dimensional complex projective space, 
class HermitianCircle extends HermitiaxiComplex2By2 : Oriented circles in 

the Riemann sphere are described by Hermitian 2x2 matrices with negative 

determinant. 

classes for numerical calculations 

package numericalMethods . function : provides interfaces which compensate the 
lack of function pointers in Java, 
interface DoubleParametrized 
interface DoubleValued 
interface DoubleArrayParametrized 
interface DoubleArrayValued 

package numericalMethods . calculus . f unctionApproximation 

class ChebyshevApproximation : provides static methods for Chebyshev in- 
terpolation and integration. 

package numericalMethods . calculus . specialFunctions 
class Clausen : implements Clausen's integral. 

package numericalMethods . calculus .minimizing : provides methods for the 
unconstrained minimization of functions of one or many variables. 

the 2D Mobius viewer 

package moebiusViewer (and sub-packages). 

Ken Perlin's 3D rcnderer 

package render : Some classes have been modified slightly by the author to allow 
resizing of the RenderApplet. The changes are documented in the code 
(search for "Springborn"). 

5.4. The class CellularSurf ace 

In appendix IdI non-oriented cellular surfaces were defined in terms of a finite 
set E of oriented edges and three permutations ((-, a, and r) of E. The class 
CellularSurf ace is an implementation of this combinatorial model. Here, we 
describe its fundamental features. See also the documentation in the javadocs 
directory. 

The oriented edges e € E are numbered consecutively, starting at 0. The edge 
numbering complies with the following convention: If edge e has number n, then 
te has number n^2 and re has number 1, where ^ denotes bitwise exclusive 
or. Hence, edge-numbers 4fc, Ak -1-1, Ak -\- 2, and 4fc -I- 3 all belong to the same 
edge of the cell decomposition. The edges with numbers Ak and Ak + 2 {Ak + 1 and 
Ak + 3) are images of each other under t, that is, they correspond to the different 
orientations of the same unoriented edge. The edges with numbers Ak and Ak + I 
{Ak + 2 and Ak + 3) are images of each other under r, that is, they correspond to 
the two preimages in the oriented double cover of the same oriented edge in the 
non-oriented cell-decomposition. 

The faces are also numbered consecutively, starting at 0, and such that if a face 
/ has number n, then the oppositely oriented face r/ has number n^l. 

The same is true for the vertices. If the vertex v has number n, then the vertex 
TV has number n^l. 
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The default constructor of CellularSurf ace produces an 'empty' cell decom- 
position with faces, vertices, and edges. 

Figure 15.21 shows the edge numbers and face numbers of the cellular surfaces 
produced by the constructors of the classes Cube and ProjectivizedCube in the 
package cellularSurf ace . examples. The first represents the combinatorial type 
of the cube. Note that orientable surfaces may be defined using only even in- 
dices. On the other side of the surface, all elements have then odd indices. The 
second represents the combinatorial type of a 'projectivizcd cube', that is, a cube 
with diametrically opposite points identified. These two CellularSurf ace s are 
constructed with calls to the method 

void buildPromFaccBoundaries(int [ ] [ ] faceBoundaries) . 

Cube is produced by a call of buildProniFaceBoundarics with the double int array 

new int [ ][ ] { 

{0, 4, 8, 12}, 
{2, 16, 32, 22}, 
{6, 20, 36, 26}, 
{10, 24, 40, 30}, 
{14, 28, 44, 18}, 
{34, 46, 42, 38} 

} 

as argument. The argument 

new int [ ][ ] { 

{0, 4, 8, 12}, 
{2, 16, 9, 22}, 
{6, 20, 13, 17} 

} 

is passed to produce ProjectivizedCube. Each row defines a face by listing the 
oriented edge numbers of the boundary in cyclic order. The first row defines face 
number 0, the second row defines face number 2, the third row defines face number 
4 and so forth. (This determines also the face boundaries of the oppositely oriented 
odd numbered faces.) The vertex indices are generated automatically. 
The methods 

int gctNumFaces() 
int getNumVertices() 
int getNumEdgcs() 

return the numbers of non-oriented faces, vertices and edges. The face indices and 
vertex indices range from to 

2 * getNumFaces() — 1 

and from to 

2 * getNumVertices() — 1 

while the edge indices range from to 

4 * getNumEdges() — 1 

because four consecutive indices correspond to a single edge of the cell decomposi- 
tion. 

Often, one wants to loop over all non-oriented faces or all non-oriented edges 
of a cell decomposition represented by a CellularSurf ace surf. This is achieved 
by loops like 
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for (int face = 0, max = 2 * surf.getNumFaces(); face < max; face 2) { 

} " 
and 

for (int edge = 0, max = 4 * surf.gctNumEdgcs(); edge < max; edge +~ 4) { 

} " 
while 

for (int edge — 0, max — A * surf.gctNumEdgcs(); edge < max; edge +— 2) { 
} " 

loops over all oriented edges of a non-oriented cell decomposition. 

Suppose an oriented edge e has index i. When called with i as argument, 

the method returns the index of (see figure ID?T|| 

int ncxtEdgcOfLcftFacc(int) the edge ere. 

int previousEdgcOfLcftFacc(int i) the edge cr~^e. 

int nextEdgcOfRightFacc(int) the edge ia'^ue. 

int prcviousEdgcOfRightFacc(int) the edge tcrte. 

int leftEdgeOflnitialVertex(int) the edge tcr~^e. 

int rightEdgeOflnitialVertex(int) the edge aie. 

int leftEdgeOfTerminalVertex(int) the edge iae. 

int rightEdgeOfTerminalVertex(int) the edge cr~^te. 

int leftFace (int) the face on the left side of e. 

int rightFace(int) the face on the right side of e. 

int initialVcrtcx (int) the initial vertex of e. 

int terminalVcrtcx(int) the terminal vertex of e. 

All these methods arc fast, because they consist essentially in an array lookup. 
There are also methods like 

int edgeWithLeftFace(int face) 

int edgeWithTerminalVertex(int vertex) 

which return an edge in the boundary of face and an edge ending in vertex^ respec- 
tively. These are computationally more expensive, because they involve a looping 
over the edges and checking, say, whether leftFace ( edge) returns face. The same 
holds for 

int boundaryLength(int face). 

To describe cell decompositions of surfaces with boundary, there is a special 
index 

static final int NO_ELEMENT = -1 

to indicate that there is no such edge/face/vertex. The class CellularSurf ace 
may be used to describe both cellular surfaces with holes and punctures and cellular 
surfaces with boundary faces and boundary vertices, as defined in section [0.21 For 
example, after calling 

buildFromFaceBoundaries(new int[ ][ ]{{0,4,8}}) 

the CellularSurf ace describes a cell decomposition of the disc with one face, three 
boundary edges, and three vertices. The methods leftFace (int) and rightFace(int) 
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The circlePattern-classes expect CellularSurf aces which encode cellular 
surfaces with boundary faces and boundary vertices. In this case every edge has a 
left face and a right face and an initial and a terminal vertex. But some faces 
may have non-closed boundary. For example, consider the example Cube (see 
figure 15.21 top) with face number 10 and its boundary edges removed. Such a 
CellularSurf ace is produced by the method call 

buildFromFaceBoundaries(new int[ ][ ] { 
{0, 4, 8, 12}, 
{2, 16, -1, 22}, 
{6, 20, -1, 26}, 
{10, 24, -1, 30}, 
{14, 28, -1, 18} 

}); 

It has 5 faces, 8 edges, and 8 vertices. The method nextEdgeOfLeftFace(int) will 
return —1 when called with argument 16, 20, 24, or 28 (and also when called with 
17, 21, 25, or 29). 
The method 

void copy(CellularSurf ace surface) 

builds a copy of surface and the method 

void buildPoincareDual(CellularSurf ace surface) 

builds the Poincare dual of surface. For example, a CellularSurf ace representing 
the combinatorial type of an icosahedron may be obtained in this way: 

import CellularSurf ace . CellularSurf ace; 
import CellularSurf ace . examples .Dodecahedron; 
CellularSurf ace surface =^ new Dodecahedron(); 
surface . buildPoincareDual( surface ) ; 

There are several methods to modify a CellularSurf ace. Figure lOl illustrates 
the simple moves. The methods 

void slideEdgeLeft(int edge) 
void slideEdgeRight(int edge) 

do not change the numbers of faces, edges and vertices, while 

int splitEdgeAlong(int edge) 

introduces a new edge and a new face, and 

int splitEdgeAcross(int edge) 

introduces a new edge and a new vertex. These last two methods return the index 
of the new edge. An edge may be contracted with 

void contractEdge(int edge). 

More complicated moves such as truncating a vertex with 
void truncateVertex(int vertex) 

and building the medial decomposition (see figure with 
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void medialGraph() 

are composed of many of the simple moves. 

5.5. The circlePattern-classes 

Here, we give a brief description of the classes which perform the construction of 
circle patterns. They are contained in the package circlePattern. variational. 
The construction proceeds in two steps. First, one of the classes EuclideanData, 
HyperbolicData, or SphericalData computes the radii of the circles, and then then 
corresponding class EuclideanLayout, HyperbolicLayout, or SphericalLayout 
computes the centers and intersection points of the circles in the respective constant 
curvature space. 

The Data classes. The classes EuclideanData. HyperbolicData, and 

SphericalData perform the minimization of the respective functionals to compute 
the radii. They all inherit from the abstract superclass GenericData which provides 
the member variables with get/ set methods. 

The protected member variables, which can be accessed with get/set methods, 

are: 

• CellularSurf ace surface. The cell decomposition which determines the com- 
binatorics of the circle pattern. Faces of the decomposition correspond to circles 
of the pattern. The surface may be non-orientable. It may be closed or a cellu- 
lar surface with boundary faces and boundary vertices (see sections 15.41 and |D.2|I . 
Boundary edges are not allowed. 

• double[ ] capitalPhi. The cone/boundary angles (f>. Before surface is set, 
capitalPhi may be null. When surface is set, capitalPhi is assigned and array with 
length surf ace. getNumFaces(). There is one angle for each non-oriented face 
/. The oriented faces with indices 2n and 2n + 1 have the cone/boundary angle 
capitalPhi [n]. 

• double [ ] theta. The exterior intersection angles 6 (see figure . Before 
surface is set, theta may be null. When surface is set, theta is assigned an array 
with length surface. getNumEdges(). There is one intersection angle 9^ for each 
non-oriented edge e. The oriented edges with indices An, An + I, An + 2 and An + 3 
have the intersection angle theta[Ti\. 

• double tolerance. The error tolerance for the minimal value of the functional. 
The default is DEFAULT_TOLERANCE. 

• double[ ] rho. The variables p. Before surface is set, rho may be null. When 
surface is set, p is assigned an array with length surface.getNumFaces(). There is 
one variable pf for each non-oriented face /. To the oriented faces with indices 2n 
and 2n + 1 belongs the variable rho[n]. Apart from the set methods, the values 
r/io[n] are changed by the method adjustRho(). 

• double value. The value of the respective circle pattern functional. This is 
set by evaluateFunctional() and adjustRho(). There is no set method. 

• double [ ] gradient. The gradient of the respective circle pattern functional. 
Before surface is set, rho may be null. When surface is set, rho is assigned an 
array with length surface. getNumFaces(). To the oriented faces with indices 2n 
and 2n + 1 belongs the gradient entry gradient[Ti\. The values gradientln] are set 
by the methods evaluateFunctional() and adjustRho(). There are no set methods. 

The functionality is provided by the following methods, which are abstract in 
GenericData: 
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• void evaluateFunctional(). From the values of surface, capitalPhi, theta, and 
rho, this method calculates the value and gradient of the functional S'cuc, Shyp, 
or iSsph (depending on the subclass of GenericData) and sets value and gradient 
accordingly. 

• void adjustRho(). Minimizes the functional S'ouc, 'S'hyp, or Ssph (depending 
on the subclass of GenericData) using a conjugate gradient method. (In the case 
of the spherical functional Ssph, the minimax procedure described in section [2. 61 is 
performed.) It attempts to determine the minimal value up to a maximal error of 
tolerance. If it returns successfully, the member variables rho, value, and gradient 
hold the approximated position of the minimum, the approximated minimal value 
and the gradient at the approximated minimum (which should be close to zero). 

• double radius(double rho). Performs the variable transformation from the 
variables p to the radii r. The relation between r and p is given by equations H2.^{|l . 
H2.14|) . or H2.23|) . depending on the subclass of GenericData. 

The Layout classes. Once the radii of the circles have been computed, 

the classes EuclideanLayout, HyperbolicLayout, and SphericalLayout compute 
the centers and intersection points of the circles. They all inherit from the abstract 
superclass GenericLayout. which provides the member variables and almost all of 
the functionality. 

The protected member variables, which can only be read with get methods, 

are: 

• ComplexProjectivel[ ] centerPoint. Complex homogeneous coordinates for 
the centers of the circles. After a call of the method layout (doia) (see below), 
centerPoint is an array of length 2*data.getNumFaces(). There is one entry for 
each oriented face. If the oriented double cover of the surface is not connected (for 
example, if the surface is orientable), then some entries of centerPoint will be null. 

• ComplexProjectivel[ ] vertexPoint. Complex homogeneous coordinates for 
the intersection points. After a call of the method layout ( rfato) , vertexPoint is an 
array of length 2*data.getNumVertices(). There is one entry for each vertex of the 
oriented double cover. If the oriented double cover is not connected (for example, 
if the surface is orientable), then some entries of vertexPoint will be null. 

• HermitiajiCircle[ ] circle. The circles of the circle pattern, described in 
a Mobius-invariant fashion as Hermitian (2 x 2)-matrices with negative deter- 
minant. After a call of the method layout ( rfaia) , circle is an array of length 
2*data.getNumFaces(). There is one entry for each oriented face. If the oriented 
double cover of the surface is not connected (for example, if the surface is ori- 
entable), then some entries of circle will be null. 

The member variables are set by the method 
void layout(GenericData data) 

of class GenericLayout. The member variables surface, theta, and rho of the 
argument data are expected to hold correct values for a circle pattern of the respec- 
tive geometry. This means, the layout method of EuclideanLayout expects to be 
passed an instance of EuclideanData with rho set by a successful call of its method 
adjustRho(), etc. 

The layout algorithm is fairly unsophisticated: Initialize member variables 
centerPoint, vertexPoint, and circle to arrays of the correct size with all entries 
equal to null. Then, place the first two points: set the centerPoint of the left 
face of the edge with index INITIAL_EDGE to the origin; and set the vertexPoint 
of the initial vertex of the same edge an appropriate distance away on the positive 
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real axis. After the initial points are placed, keep searching for more points whose 
position can be determined until no more can be found. 

Almost all the work is done in the abstract superclass GenericLayout. Its 
layout method is not abstract. But internally, it calls the two abstract methods of 
GenericLayout: 

abstract void assignOthcrFixPoint(ComplexProjectivel center, 

ComplexProjectivel otherFixPoint); 
abstract void assignCirclc(HermitiaiiCircle circle, 

ComplexProjectivel center, 

double radius); 

All the subclasses EuclideanLayout, HyperbolicLayout, and SphericalLayout 
do is to implement these two abstract methods: 

The method assignOtherFixPoint assigns the argument otherFixPoint depend- 
ing on the argument center. In the class EuclideanFunctional, the argument 
otherFixPoint is always set to the infinite point of the euclidean plane. In the 
class HyperbolicFunctional, otherFixPoint is set to image of center under an 
inversion on the circle which is the boundary of hyperbolic space. In the class 
SphericalFunctional, otherFixPoint is set to the point which is diametrically 
opposite to center on the Riemann sphere. 

The method assignCircle assigns circle by the circle which, in the respective 
geometry, has center center and radius radius. 

5.6. Computing Clausen's integral 



In 1832, Clausen first tabulated the integral l|B.3(l . which now bears his name 
|15| . He expands the functions 

Cl2(7ra;) + ^(2 log + xlog ((2 + x)(2 - x)\x\)) 

and 

1 + X 

Cl2(7r(l - x)) - TTlog TTX-log ((1 + a;)(l - x)) 

1 — X 

into power series of x, which he uses to calculate Cl2(x) in the intervals [0, ^] and 
[|,7r] to 16 decimals. 

We proceed in a similar fashion, except that we use Chebyshev series. For 
X g [— 7r,7r], the function 

, , /2sin§' 

h{x) = - log 

V 2; 

is analytic and 

Cl2(x)= r h{Odi^x{\og\x\~\). 







The function h{x) is approximated in the interval [— tt, tt] by a Chebyshev series, 
which can easily be integrated. 

The algorithms for fitting, integrating and evaluating Chebyshev series are 
described in the Numerical Recipes |33j . They are implemented in static meth- 
ods of the class ChebyshevApproximation, which is contained in the package 
numericalMethods . calculus . f unctionApproximation. 

Using these methods for Chebyshev approximation, the class Clausen in pack- 
age numericalMethods . calculus . specialFunctions provides the static method 

double cl2(double x) 

to evaluate Clausen's integral. It also provides Catalan's constant G = Cl2(f ): 
static final double CATALAN 
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Figure 5.2. Cube and ProjectivizedCube. 
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Figure 5.3. Simple moves. 



APPENDIX A 



Proof of the trigonometric relations of lemma 12.51 

and lemma 12.101 



In this appendix, we derive the below formulae for the remaining angles in 
a hyperbolic or spherical triangle when two sides and the enclosed angle between 
them are given. Suppose the given sides of the triangle are ri and r2 , the included 
angle between them is 6, its complement is 6* = ti — 9, and the remaining angles 
are ipi and (/32, as shown in figure 

In the spherical case, 



1 1 + eP^^P^'^ 



— log 

2i 1 - ePi+P^-^e* 



(A.l) 



where we assume < r , < tt, and 



P3 



log tan . 
^ 2 



In the hyperbolic case. 



"Pi 



1 , 1 
— log — 
2i ^ 1 



, gP2-pi+i9' \ _j_ gPl+P2+if' 

— — log 

, gP2-pi-*f* 2i 1 + ePi+P2-*''* 



(A.2) 



where 



P3 



log tanh — 
^ 2 



These equations (|A.1|I and (|A.2(I imply equations (|2.22() and (|2.13|1 . 

To derive such formulae from known relations of non-euclidean trigonometry 
can be cumbersome. Instead, we present a method using Mobius transformations 
which has proved to be helpful. We identify the group of Mobius transformations 
of C = C U {oo} with the projectivized linear group PGL{2, C). A Mobius trans- 
formation 



az + b 
cz + d ^ 



corresponds to 



:= C 



ad-bc^O 



e PGL{2,C). 



A.l. The spherical case 

The unit sphere = {{xi,X2,X3) € M.^ | xi^ + + x^'^ = l} is mapped to 
the extended complex plane C = C U {oo} by stcrcographic projection 

Xi+ 1X2 

[Xi,X2,X3) ^ Z = — . 

1- X3 

The orientation preserving isometrics of the sphere (the rotations) correspond to 
Mobius transformations ^] with a = d and b = — c. Rotations around the xi-, 
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X2-, and X3-axis by an angle a correspond to 
R2{a) = 



— sm 
cos ^ 



g-ia/2 



i?2(a) 



I tan ^ 



tan ^ 





e*" 0" 







1 


then 






i tan 


Q 

2 




1 




7 


- tan 


Q " 

2 




1 







The quantities ri, r2, I, and (^i, (^2, ^ = tt — 9* are the sides and angles of a spherical 
triangle as in figure if and only if 

i?2(r-i)i?3(e*)i?2(r-2) = i?3(^l)i?2(0^3(<^2). 
With = logtan(rj/2), the left hand side of this equation equals 

" 1 -eP^ 
ef^ 1 



LHS 









■ 1 -eP^' 




1 




eP^ 1 



gPl+P2 



oP2 



0PI+P2 



_gP2+iS _ gPl 

_ gPl+P2+i9* 



„P2+ie 



+ e 



Pi 



1 



oPl+P2+i£' 



1 



„Pl+P2+iS 



■ e 



P2 



1 



„pi+P2+4e 



With A = logtan(^/2). the right hand side equals 
RHS = 



Hence, 







■ 1 -e^' 




'e^^^ 0" 


1 




1 




1 










1 





„P2+j6(* 



oPi 



cpi = arg ■ 



oPi+P2+ie 



which implies equation (|A.ip . 

A. 2. The hyperbolic case 

We use the Poincare disc model of hyperbolic 2-space: = £ C \ \z\ < 1^ 
with metric \ds\ = 2\dz\/{l — The orientation preserving isometrics of are 

the Mobius transformations [ ° ^ ] with a = d and b = c. Let 

3sh I sinh | _ 1 tanh : 
nh I cosh I tanh | 



R{a) 



1 



Ja/2 










e*" 0" 


g-iQ/2 




1 



The transformation T(s) translates all points on the line £ H^\lmz = 0} 
towards z = 1 by the distance s. The transformation R{a) is a rotation around 
2; = by the angle a. 
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The quantities ri, r2, I, and ipi, ip2, 9 = tt — 9* are the sides and angles of a 
hyperboUc triangle as in figure if and only if 

T{ri)R{9*)T{r2) = R{ipi)T{l)R{ip2). 

With pj = logtanh(rj /2), the left hand side of this equation equals 



LHS 



1 







1 

= P2 



oP2 



0PI+P2 



0PI+P2 



1 

eP2+ir ^ 
gP2+»r 



,pi 



1 



oPl- 



Pl+P2+'t9* 

„eP2 



1 



1 



Pl+P2+iS' 



Pl+p-2+i6' 
1 



With A = logtanh(//2), the right hand side equals 
RHS = 



Hence, 







"1 e^" 
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oPi 



If I = arg ■ 



1 + e 



Pi+P2+ie* ' 



which implies equation 



APPENDIX B 



The dilogarithm function and Clausen's integral 



In this appendix, wc collect all relevant facts about the dilogarithm function 
and Clausen's integral that is relevant for this paper. A more thorough treatment 
and an extensive bibliography are contained in Lewin's monograph |29j . 

For \z\ < 1, the dilogarithm function is defined by the power series 

^3 



For \z\ < 1, 



and hence 



Ll2(^) = Y^ + ^ + ^ 



z z z 
log(l-.) = - + y + y 



Li.W^-fi^^fc^dC. (B.l) 
Jo s, 

In the light of this integral representation, one sees that the dilogarithm can be 
continued analytically to the complex plane cut from 1 to c» along the positive real 
axis. 

Take logarithms on both sides of the identity 

1 - z" = (1 - z)(l - ujz){\ - u?z) ... (1 - cj"-^z), 

where a; = e^^*/" is the fundamental n*'^ root of unity, then divide by z and integrate 
to obtain 

- Li2(z") = Li2(z) + Li2(tJz) + Li2(tJ^z) + . . . Li2(a'"~^z) (B.2) 
n 

for l^l < 1. Both sides of the equation can be continued analytically to the complex 
plane with radial cuts outward from the n**^ roots of unity to infinity. 

Clausen's integral Cl2(a;) can be defined by the imaginary part of the diloga- 
rithm on the unit circle: 

Cl2(x) = ImLi2(e") 

= i(Li2(e-)-Li2(e-^-)). 

(The name "Clausen's integraV comes from the integral representation (IB.BI) .) We 
consider Clausen's integral as a real valued function of a real variable. It is 27r- 
periodic and odd. The power series representation of the dilogarithm yields the 
Fourier series representation for Clausen's integral, 

ci2(x)=f:^. 

n=l 

Substitute ^ = e*^ in the integral representation of the dilogarithm (|B.ip to obtain, 
for < X- < 27r, 

Cl2(x) = - riog('2sin0 dC. (B.3) 
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By periodicity, for all x e R, 

^Cl2(x) = -log|2sinf |. (B.4) 
Clausen's integral is almost the same as Milnor's Lobachevski function |31| 



Jl{x) 



^"log|2sin^|de= ici2(2x) 



Lobachevski |30| himself defined a function L{x) (Lobachevski's Lobachevski func- 
tion) by 



Lix) = - / log(cosOd^, 
Jo 



such that 

Lix) = -Cl2(2a;-7r) + a;log2. 



From equation l|B.2|) one obtains 

n — 1 

- CUnx) = V Cl2(a; + 27rfc/n) . 

77, ^ — ^ 



fe=0 

In particular, using that Cl2(a;) is 27r-periodic and odd, one obtains the double-angle 
formula for Clausen's integral: 

ici2(2a;) = Cl2(a;)-Cl2(7r-a;). (B.5) 

We will now derive a formula expressing the imaginary part of the dilogarithm 
in terms of Clausen's integral not only on the unit circle, but anywhere in the 
complex plane. Suppose x is real and < 9 < 2tt. Substitute C = e^+*^ in (|B.ip . 
to obtain 

ImLi2(e-+'^'') = l(Li2(e^+''') - Li2(e^"*«)) 

1 r /l-e^~*^\ 

Now substitute 

1 , /I - e^-'^ 



(B.6) 



2i V 1 - 
and note that inversely 

2 sin?7 



2 sin(7/ + 6) 
Finally, one obtains 

ImLi2(e^+^«) ^yx+^ Ch{2y) - ^ Ch{2y + 29) + ^ Ch{20) , (B.7) 

where 

From this, we derive the formula 
ImLi2(e^+*'') +ImLi2(e--^+*^) = 

px + ci2(p + r ) + ci2(-p -I- e*) - ch{2e*) , (b.8 

where 6* = ir — 6, and 



1 (1 + 6^+^*' ){l + e-^^'^ ) 
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Finally, p and x are related by 



tan I ^ ) = tanh ( — ) tan ( — 



APPENDIX C 



The volume of a triply orthogonal hyperbolic 
tetrahedron with a vertex at infinity 

Milnor |31l pp. 19f| calculated the volume of a triply orthogonal hyperbolic 
tetrahedron (three-dimensional hyperbolic orthoscheme) with two vertices at in- 
finity by a straightforward integration. By the same computation, one derives a 
formula for the volume of a triply orthogonal hyperbolic tetrahedron with one ver- 
tex at infinity. For reference, we present this computation here. 

Figure [T^. II shows a triply orthogonal tetrahedron with one vertex at infinity in 
the Poincare half-space model of hyperbolic space. In the figure, the tetrahedron is 
truncated by a horosphere centered at the infinite vertex. Three of the six dihedral 
angles are ^. Because the sum of dihedral angles is tt at an infinite vertex and 
larger than tt at a finite vertex, the three remaining angles are a, ^ — a, and /3, as 
shown in the figure, where 

< a < f3 < -. 
- 2 

If a = /?, the tetrahedron has two infinite vertices. 

Proposition. The hyperbolic volume of the triply orthogonal hyperbolic tetra- 
hedron with a vertex at infinity is 

y = i (2 Cl2(7r - 2a) + Cl2(2a - 2/3) + Cl2(2a + 2/3)) . (C.l) 

Setting a = /3 and using the double-angle formula for Clausen's integral (|B.5|l . 
one obtains the formula for the volume of a triply orthogonal hyperbolic tetrahedron 




Figure C.l. Left: A triply orthogonal hyperbolic tetrahedron with one vertex at 
infinity (truncated at a horosphere). Right: The domain of integration in the xy- 
plane. 
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with two vertices at infinity and characteristic angle a: 



V 



ici2(2«). 



(C.2) 



Proof of the proposition. In the Poincare half-space model, hyperbolic 
space is = {{x, y, z) £ M."^ \ z > 0} with metric ds"^ = -^{dx'^ + dy'^ + dz"^). The 
volume form is therefore \dx dy dz. Hence, 



V 



cos /3 



x=0 



V=0 



„2 Z 



— dz dy dx, 



see figure (right). The first two integrations are readily performed: 



V = - 



cos /3 



2 ■'x=o 



1 



■ dy dx 



cos /3 



V=0 
1 



1 — x'^ — y 

( xtana 



Substitute x 



lx=0 

cost, and note that 



: artanh 



dx. 



artanh 



tana 
tant 



1 2.-1- '^^^ ^ cos t 
Yqo- COS a sin t 

2 \ sin a cos t 

cos a sin t 

COS a sin t + sin a sin t 



Using equation (jB.Sp . one obtains 



V 



1 
1 r 



t: log 

2 cos asmt- 

1 , sin(t + a) 
— log — ^ 

2 ^ sin(t - a) 



sm(t + a) 

log ^77 7 

sm(t — a) 



sin a sin t 



Cl2(2t-t-2a)-Cl2(2t-2a) 



and hence equation IjCljl . 



APPENDIX D 



The combinatorial topology and homology of 
cellular surfaces 

This appendix is concerned with a combinatorial model for finite cell decompo- 
sitions of surfaces, and with their Z2-homology. Equivalent combinatorial models 
were described by Jacques (23, , Tutte j46| , and Jackson and Visentin |22) . They are 
sometimes called 'winged-edge models'. We present as much of the Z2-homology 
theory as needed in the proof of lemma [^.141 in section It). 3. 51 Of course, all of this 
is well known. 



Consider cell decompositions of a compact oriented 2-manifold without bound- 
ary. There is one such cell decomposition with no edges: The sphere with one face 
and one vertex. The others can be described by two permutations of the oriented 
edges: The first, t, reverses the orientation of an edge. The second, cr, maps each 
oriented edge to the next oriented edge in the boundary of the face on its left side. 
(See figure EU) 

Definition. An oriented cellular surface S = {E, l, a) is a finite set E together 
with an involution l : E —y E without fixed points (that is, = id and lc ^ e for 
all e G E) and a permutation a : E ^ E. 

The elements e £ E arc the oriented edges of E. The unoriented edges are 
unordered pairs |e | = {e, te}. 

The faces of S are the orbits in E of the group generated by a. The face on 
the left side of an oriented edge e £ E is the orbit of e. The face on the right side 
of e is the orbit of be. 

The vertices of E are the orbits in E of the group generated by l(t~^. The 
initial vertex of an oriented edge e G -E is the orbit of e. The terminal vertex of e 
is the orbit of te. 

An oriented cellular surface is connected if the group generated by l and a acts 
transitively on e. 



Figure D . 1 . Elements of the group generated by the permutations l and t acting 
on an oriented edge e. 



D.l. Cellular surfaces 
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Figure D.2. How the 'winged edges' are glued together. 

Every oriented cellular surface corresponds to a cell decomposition of a com- 
pact oriented 2-manifold without boundary. Indeed, the cell decomposition can be 
constructed from the cellular surface as follows. For each oriented edge e G E, take 
one triangle (a 'winged edge') and glue them together as shown in figure lD?2l 

Definition. Let Ei = {Ei,Li,ai) and E2 = (£^2, ^2, ^2) be two oriented cel- 
lular surfaces. An orientation preserving homeomorphism f from Ei to E2 is a 
bijection f : Ei ^ E2 with 

f O Ll = L2 O f and / O O"! = (72 ° /• 

An orientation reversing homeomorphism f from Ei to E2 is a bijection f : Ei E2 
with 

/oti=i2o/ and / o (lO-f ^ti = (T2 o /. 

An orientation preserving homeomorphism / maps the orbits of ii, cti, and 
ii o to the orbits of L2, (72, and t2 ° <^2^ ■ Hence, it induces bijections between 
the non-oriented edges, the faces, and the vertices of Ei and E2. The face on the 
left side of e G Ei is mapped to the face on the left side of /(e), and the initial 
vertex of e is mapped to the initial vertex of /(e). 

An orientation reversing homeomorphism / also induces bijections between the 
non-oriented edges, the faces, and the vertices of the two cellular surfaces. The 
orbits of Ll and ii o crj"^ are mapped to the orbits of i2 and t2 ° (although the 
latter are traversed in the opposite direction). The orbits of ui are not mapped to 
orbits of 0-2, but to the images of such orbits under L2- Nonetheless, this puts the 
faces in a one-to-one correspondence: The orientation reversing homeomorphism / 
maps the face on the left side of e € to the face on the right side of /(e). 

Cell decompositions of non-oriented (and possibly non-orientablc) 2-manifolds 
are described in terms of their orientable double cover. 

Definition. A non-oriented cellular surface E — [E, t, c, r) is an oriented cel- 
lular surface Eq = (-B, t, a) together with an orientation reversing homeomorphism 
r from Eg onto itself, which is an involution without fixed oriented edges. The faces, 
oriented edges, and vertices of E are pairs of faces, oriented edges, and vertices of 
Eq, which are mapped to each other by r. 

In other words, t : E E is another permutation with = id, T(e) 7^ e for 
all e^E, TL~ LT, and na^^L = ar. 

Definition. The Poincare dual cellular surface of an oriented cellular surface 
E — {E,L,a) is the oriented cellular surface E* = {E,L,a*) with a* = La~^. The 
Poincare dual cellular surface of a non-oriented cellular surface E = {E,L,a,T) is 
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the non-oriented cellular surface S* = {E, l, a* ,t*) with a* as in the oriented case 
and T* = LT. 

There is a one-to-one correspondence between the unoriented edges of E and 
E*. between the vertices of E and the faces of E*. and between the the faces of E 
and the vertices of E*. The Poincare dual of the Poincarc dual is essentially the 
same cellular surface. The map e'—^ be is an orientation preserving homeomorphism 
from E to E**. 

D.2. Surfaces with boundary 

Combinatorial models for cell decompositions of surfaces with boundary are 
less canonical. In this section we describe a few ways in which the above definition 
for cellular surfaces may be adapted to deal with bounded surfaces. We treat only 
the oriented case; unoriented surfaces may be described in terms of their oriented 
double cover, as in the previous section. 

The most straightforward way to define oriented cellular surfaces with boundary 
is to mark some faces as holes: 

Definition. An oriented cellular surface with holes E = {E, l, a, H) is an 
oriented cellular surface {E, l, a) together with a subset H of the set of faces, such 
that, if an oriented edge e € E is contained in a face of H, then te is not contained 
in a face of H. 

There is no Poincarc duality for cellular surfaces with boundary by this defini- 
tion. As a remedy, one may consider surfaces with holes and punctures: 

Definition. An oriented cellular surface with holes and punctures 

j:^{E,L,a,H,P) 

is an oriented cellular surface {E, l, a) together with a subset H of the set of faces 
and a subset P of the set of vertices, such that, if an oriented edge e G E is contained 
in a face of i7 or a vertex of P, then te is not contained in a face of H or a vertex 
in P. 

(If Poincare duality is to be upheld, and adjacent holes are not allowed, then 
adjacent punctures must be forbidden, too. Alternatively, one may allow adjacent 
holes and punctures.) 

In the computer implementation described in chapter a slightly more exotic 
combinatorial model is used. It is obtained by deleting all boundary edges. As a 
consequence, some faces have non-closed boundary; and therefore, for some oriented 
edges, there is no next edge in the boundary of the left face. This model is well 
suited for applications where one is not interested in boundary edges at all, as in 
the case of circle patterns. 

Definition. An oriented cellular surface with boundary faces and boundary 
vertices E = {E, l, a) is a finite set E together with an involution l : E ^ E 
without fixed points and a map a : E ^ E U {0} which satisfies the following 
condition: 

V 6*1, (3*2 e -E : cr(ei) = cr((3*2) 7^ ei = 62. 

Although a is not invertible, there is a 'pseudo-inverse' <t^^ : E ^ E U {0} 
which satisfies a^^{a{e)) = (3* for aile £ E with a{e) ^ 0, and aiu^^ie')') = (3* for all 
cgE with a-~^{e) ^ 0. Faces and vertices are defined as 'orbits' of cr and tcr~^ as 
in the case of ordinary oriented cellular surfaces. Every edge has a left and a right 
face and an initial and a terminal vertex, and the Poincare dual cellular surface is 
well defined. 
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D.3. Z2-Homology 

Let E be an oriented cellular surface, or, equivalently, a cell decomposition of 
a compact 2-manifold without boundary. Let F, E, and V be the sets of faces, 
non-oriented edges, and vertices. 

D.3.1. Homology groups. The chain-spaces Co(S), Ci(S]), and C2(S) are 
the 1i2 vector spaces generated by V, E, and F; that is, they arc the vector spaces of 
formal sums of vertices, edges, and faces with coefficients in the field I12 — TLjTL. 
The boundary operators d : C2(E) — > Ci{T,) and d : Ci{T.) — > Co(E) are the 
Z2-liiiear maps which are defined as follows. A face / G F of a cellular surface 
was defined as an orbit / = {ei, . . . , e„} of oriented edges. The boundary of / is 
9/ = |ei| + . . . + |e„|, where |e | is the non-oriented edge corresponding to e. For 
an edge e £ E with end vertices vi,V2 G V , de = vi-\- V2- 

The chain complex Ch.(S]) is the following sequence of chain spaces and linear 
maps: 

^ C2(S) ^ ^ Co(E) ^ 

The cycle spaces are the kernels Z„{Y,) = {c e C„(S) | 9c = 0} for n £ {0, 1, 2}. The 
boundary spaces are the images i?2(E) ~ and i3„(E) = (9C„+i(I]) for n € {0, 1}. 
The homology groups are the quotient spaces Hn(T,) = Z„(S])/i?„(I]) and their 
dimensions are the Betti numbers /i„(S) = dim^a -ff)i(E). 

The second Betti number /12(E) is the number of equivalence classes of faces 
under the equivalence relation 

. . Every closed 2-chain which 

contains /i also contains f2- 

The zeroth Betti number /io(E) is the number of equivalence classes of vertices 
under the equivalence relation 

The vertices z'l and V2 can be 
jomed by a path. 

Using subdivision arguments one can show that /12(E) = /io(E). Of course, this is 
just the number of connected components of E. 

D.3.2. Poincare duality. The Poincarc dual cellular surface E* has face set 
F* = V, edge set E* = E, and vertex set V* ~ F. The chain spaces are therefore 
C„(E*) = C2-n(E) for n G {0,1,2}. Denote the boundary operators of the chain 
complex C*(E*) by d* . Let the chain spaces be equipped with scalar products 
( • , • ) such that F, E, and V are orthonormal bases. 

Lemma D.l (Poincarc duality). The boundary operators d* of the chain com- 
plex C*(E*) are the adjoint operators of the boundary operators d of the chain 
complex C* (E) with respect to the scalar product ( • , • ) . 

This means, for f e F = V* , e e E ^ E* , and v e V ^ F* , 

{df,e) = {f,d*e) and {de,v) = {e,d*v). 

Proposition (Poincarc duality for cellular surfaces). The Betti numbers ofT, 
and 12* satisfy 

/i„(E*) = /i2-„(E). 

Proof. The dimensions of i/„(E*) and i/2-n(E) are equal because they are 
dual vector spaces. This follows from the following facts from linear algebra. 

Let y be a vector space, and let U and W be subspaces such that W C U C V. 
Then the dual space of the quotient U /W is 

{u/wy = w^/u^, 
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C2(T) Ci(T) — ^ Co(T) 



> C2(E) Co(S) 



C2(S,T) Ci(S,T) — ^ Co(S,T) . 















Figure D.3. The short exact sequence -> C,(T) C,(S) 
chain complexes. 



C,(E,T) ^ of 



where C V* and C/^ C V^* arc the spaces of hncar functionals on V which 
vanish on W and U, respectively. 

Let / : V — > be a linear map between two vector spaces, and let /* : 
W* — > V* be the adjoint map. Then the kernels and images of / and /* are 
related by 

Imr = (Ker/)^, Kerr = (Im/)^. 

Hence, 

i7„(S)* = (Z„(E)/B„(E))* ^ B„(S)^/Z„(S)^ = 

= Z2-„(E*)/S2-„(S*) = JJ2-„(S*). 



D.3. 3. Relative homology. Define a subcomplex T of a cellular surface E as 
a triple T = (Ft, St, Vt) of subsets Ft C F, Ft C F, and Vr C V, such that the 
following holds: If a face in Ft contains an oriented edge e, then the corresponding 
unoriented edge |e | £ Ft; and the end- vertices of any edge in Ft are contained in 
Vr- 

The relative chain complex C*(E, T) consists of the quotient spaces C„(E, T) = 
C„(E)/C„(T) and the boundary operators induced by the boundary operators of 
C*(S). The corresponding homology groups _ff„(E, T) arc the relative homology 
groups of the pair (E, T). 

The composition of the inclusion map C„(T) C„(E), followed by the natural 
projection C„(E) — > C„(E,T) is the zero map. Hence, the short sequence 

— > C,(T) a(E) — > C,(E, T) — >0 

is exact. This means that the columns of the commutative diagram shown in fig- 
ure lD.3l are exact. By the zigzag-lemma, there is a long exact sequence of homology 
groups 

^ F2(T) F2(E) F2(E, T) 

-^i/i(T) ^Fi(E) ^iJi(E,T) 

Ho{T) Fo(E) iJo(E,T) 0. 
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— ^ Ci(r) — ^ Co(r) > 



C2(S) Co(S) . 



^ C2(s,r) — ^ Ci(E,r) Co(s,r) . o 





Figure D.4. The short exact sequence ^ C^T) C*(S) ^ C*(E, T) ^ 0. 



D.3.4. Lefschetz duality. Let T = (Fx,i?T,^T) be a subcomplex of the 
cellular surface S as in the previous section. Then the complement of T in S*, 

S* \ T {F* \ Vr, E* \ Er, V* \ Fr), 

is a subcomplex of S*. 

Proposition (Lefschetz duality for cellular surfaces). The Betti numbers of 
the subcomplex S* \ T and the pair (S, T) satisfy 

/i„(I]*\T) = /i2-„(S,T). 

Proof. As in the case of Poincare duality, the dimensions of the homology 
groups \ T) and i?2-ri(S, T) are equal because they are dual vector spaces. 

This is proved in the same way. ■ 



D.3.5. Proof of lemma f2.14[ Consider the long exact homology sequence of 
relative homology for a 1-dimcnsional subcomplex P, that is, a cellularly embedded 
graph. Let Er C E and C F be the edge and vertex sets of P. Since C2(P) is 
zero, the short exact sequence 











a(P) ^C7,(S) ^C,(S,P) 

of chain-complexes is equivalent to the commutative diagram shown in figure ID^ 
in which the columns are exact. This gives rise to the following long exact homology 
sequence. 

F2(S) H2{^,r) Fi(P) J?i(E) ffi(S,P) 

^ HoiT) Hoi^) i/o(S,P) 
Since the sequence is exact, the corresponding dimensions satisfy 

/i2(s) - /i2(s, p) + /ii(p) - + p) - ho{r) + ho{j:) - /io(s, p) = o. 

Now assume that S is connected and P is non-empty, that is, Vr 7^ 0- Then 
/io(S,P) = 0, because from every vertex of S there is a path to a vertex of the 
graph P. Also, 

hi{T)~ho{r) = \Er\~\Vr\ 

and 

h2{^) - + /^o(S) = l^^l - \E\ + \Vl 
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and hence, 

/i2(S,r) - \Er\ + \Vr\ = |F| - \E\ + \V\ + /ii(S,r). 
Now, /i2(S, r) is the number of regions into which the graph F separates the cellular 
surface S. By Lefschetz-duality /ii(I],r) = \ F). Suppose T separates S* 

into the components IIi, . . . , n„; then \ T) = ^ /ii(nj). 

This completes the proof of lemma 12. 141 
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